Mon. Not. R. Astron. Soc. 000,|TH23l(2010) Printed 14 June 2010 (MN ET^ style file v2.2) 



O 

3 



The LABOCA survey of the Extended Chandra Deep Field South: 
A photometric redshift survey of submillimetre galaxies 

J. L. Wardlow 1 *, Ian Smail 2 , K. E. K. Coppin 2 , D. M. Alexander 1 , W. N. Brandt 3 , A. L. R. 
Danielson 1 , B. Luo 3 , A. M. Swinbank 2 , F. Walter 4 , A. WeiB 5 , Y. Q. Xue 3 , S. Zibetti 4 , F. 
Bertoldi 6 , A. D. Biggs 7 , S. C. Chapman 8 , H. Dannerbauer 4 , J. S. Dunlop 9 , E. Gawiser 10 , 
R. J. Ivison 9 ' 11 , K. K. Knudsen 6 , A. Kovacs 5 , C. G. Lacey 2 , K. M. Menten 5 , N. Padilla 12 , 
H.-W. Rix 4 , and R R van der Werf 13 



^Department of Physics, Durham University, South Road, Durham, DH1 3LE, UK 
2 Institute for Computational Cosmology, Durham University, South Road, Durham, DH1 3LE, UK 

3 Department of Astronomy and Astrophysics, 525 Davey Lab, Pennsylvania State University, University Park, PA 16802, USA 
^Max-Planck-Institut fur Astronomie, Konigstuhl 17, D-691 17 Heidelberg, Germany 
6 ' Max-Planck-lnstitut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 
6 Argelander lnstitut fiir Astronomic, Auf dem Hiigel 71, D-53121 Bonn, Germany 
7 European Southern Observatory, Karl-Schwarzschild Strasse 2, D-85748 Garching, Germany 
^ ^) i ^Institute of Astronomy, Madingley Road, Cambridge, CB3 OHA, UK 

• 9 SUPA (Scottish University Physics Alliance), Institute for Astronomy, University of Edinburgh, Royal Observatory, Edinburgh EH9 3HJ, UK 
■ 10 Physics and Astronomy Department, Rutgers University, Piscataway, NJ 08854, USA 
Oh- 11 UK Astronomy Technology Centre, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 

12 Departmento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, Santiago, Chile 
^Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 



o 



I 

o 



CO 



14 June 2010 



' ABSTRACT 



(N 

o 
o 

> 



We derive photometric redshifts from 1 7-band optical to mid-infrared photometry of 74 robust 
radio, 24 |xm and Spitzer IRAC counterparts to 68 of the 126 submillimetre galaxies (SMGs) 
selected at 870 |im by LABOCA observations in the Extended Chandra Deep Field South 
(ECDFS). We test the photometric redshifts of the SMGs against the extensive archival spec- 
troscopy in the ECDFS. The median photometric redshift of identified SMGs is z = 2.2 ±0.1, 
the interquartile range is z = 1.8-2.7 and we identify 10 (~ 15%) high-redshift (z ^ 3) 
SMGs. We derive a simple redshift estimator for SMGs based on the IRAC 3.6 and 8 \xm 
fluxes which is accurate to Az ~ 0.4 for SMGs at z < 4. A statistical analysis of sources 
around unidentified SMGs identifies a population of likely counterparts with a redshift distri- 
bution peaking at z = 2.5 ± 0.3, which likely comprises ~ 60% of the unidentified SMGs. 
This confirms that the bulk of the undetected SMGs are co-eval with those detected in the 
■ radio/mid-infrared. We conclude that at most ~ 15% of all the SMGs are below the flux lim- 

its of our IRAC observations and lie at z > 3 and hence around ~ 30% of all SMGs have 
z > 3. We estimate that the full SWo^m > 4 mJy SMG population has a median redshift of 
2.5 ± 0.6. In contrast to previous suggestions we find no significant correlation between sub- 
millimetre flux and redshift. The median stellar mass of the SMGs derived from SED fitting 
is (9.2 ± 0.9) x 10 10 Mq and the interquartile range is (4.7-14) x 10 10 Mq, although we cau- 
tion that the uncertainty in the star-formation histories results in a factor of ~ 5 uncertainty 
in these stellar masses. Using a single temperature modified blackbody fit with j3 = 1.5 the 
median characteristic dust temperature of SMGs is 35.9 ± 1.4K and the interquartile range 
is 28.5^-3.3 K. The infrared luminosity function shows that SMGs at z = 2-3 typically have 
higher far-infrared luminosities and luminosity density than those at z = 1-2. This is mir- 
rored in the evolution of the star-formation rate density (SFRD) for SMGs which peaks at 
z — 2. The maximum contribution of bright SMGs to the global SFRD (~ 5% for SMGs with 
SWo^m ~ 4mJy; ~ 50% for SMGs with Ss7o^ m > 1 mJy) also occurs at z ~ 2. 
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1 INTRODUCTION 

Observations in the millimetre and submillimetre wavebands pro- 
vide a uniquely powerful route to survey the distant Un iverse for in- 
tense dust-obscured starbursts dB lain & Longairll 1993b . This is due 
to the negative K-correction arising from the shape of the spectral 
energy distribution (SED) of the dust emission in the rest-frame 
far-infrared, which results in an almost constant apparent flux for 
sources with a fixed luminosity at z ~ 1-8. 

Over the past decade, a series of ever larger surveys in the 
submillimetre and millimetre wavebands have mapped out a popu- 
lation of sour ces at mjy-fl ux limi ts with a surprising l y high surface 
densi t y (e.g. ISmail et alJ 1 19971; iB arger et al. 1998; Hughes et al. 
19981 : lEales et all 1 19991: iBertoldi et al.1 |200CL Ho07l: ICop pin etafl 
20061: [Knudsen et al. 2008; IWeifi et alJ 120091 ; [Austerman n et al.1 



2010). The mjy fluxes of these sources imply far-infrared lumi- 
nosities of > 10 12 L©, if the sources are at cosmological distances, 
z > 1, classing them a s ultraluminous infrared galaxies (ULIRGs; 
I Sanders & Mirabellll996f) . Their high surface density is far in ex- 
cess of that expected from a "no evolution" mod el, suggesting very 
strong evolution of the population: oc (1 + z) 4 (Smail et al. 1997; 
iBlainetalJ 1999b . If this results from strong lumin osity evolution of 
starburst galaxies (as opposed to obscured AGN; Alexand er et al.l 
120051 then a significant fraction of the massive star formation (and 
metal production) at high redshift may be occurring in this popula- 
tion. 

To confirm this evolution and understand the physical pro- 
cesses driving it requires redshifts for the submillimetre galaxies 
(SMGs). Due to the coarse spatial resolution of the submillimetre 
and millimetre maps from which the SMGs can be identified, com- 
bined with their optical faintness (in part due to their high dust ob- 
scuration), it has proved challenging to measure their spectroscopic 
redsh ift distribution (e.g. lBarger et al.ll999l : IChapman et ail 2003a. 
120051) . 

In fact, spectroscopic redshifts are not necessary to map 
the broad evolution of the SMG population and cruder photo- 
metric redshifts can be sufficient, if they are shown to be re- 
liable. Various photometric redshift techniques have therefore 
been applied in an attempt to trace the evolution of SMGs, 
using their optical/near-/mid-infrared or far- i nfrared/radio SEDs 
(e.g. ICarilli & Yurj [l999l: ISmail et al.1 l2000l; llvison et all 1 2004 ; 
Pope et all 12001 1 20061: llvison et al. 20071: lAretxaga et al] 12007 ; 



Clements et al. 2008: 



Dye et alfeOollBiggs et alj|2010h 



Both spectroscopic and photometric analyses suggest that the 
bulk of the SMG population lies at z > 1, with an apparent peak 
at z ~ 2.2 for the subs et of SMGs which can be located through 
their /ijy radio emission dChapman et alj2005h . Nevertheless, there 
are significant disagreements between the different studies (se e 
e.g. IChapmanetaflfeoOSl : Iciements etafl 120081 : iDve et al J 120081) . 
which may arise in part due to differing levels and types of in- 
completeness in the identifications and biases in the redshift mea- 
surements. The most serious of these is the incompleteness due to 
challenges in reliably locating the correct SMG counterpart. They 
are typically identified through statistical arguments and physical 
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these locate only ~ 60-80% of SMGs. The expectation is that the 
SMGs whose counterparts are missed could potentially include the 
highest redshift (and thus the faintest in the radio and mid-inf rared) 
examples, biasing the derived evolution dlvisonetalj|2005l) . At- 
tempts to address this incompleteness through time-intensive sub- 



millimetre interferometry have l ocated a small fraction o f pre- 
viously unidentified SMG s (e.g. iDannerbauer et all I2002L 120081 : 
lYounger etail 12001 120091 ; IWang et alj|2007l) but the nature and 
redshifts of this unidentified subset of SMGs remains a critical is- 
sue for studies of the population as a whole. 

In this paper we use optical, near- and mid-infrared pho- 
tometry to study SMGs detected in the Extended Chandra Deep 
Field Sout h (ECDFS) by the Large APEX BOlometer CAmera 
(LABOCA: lsiringo et : al.l2009h o n the Atacama Pathfinder Exper- 
iment (APEX; lGUsten et al.ll2 006) 12-m telescope in the LA BOCA 
ECDFS Submillimetre Survey (LESS; IWeifi et al.1 l2009h . LESS 
mapped the full 30' x 30' ECDFS at 870- iim to a noise level of 
o"870(im ~ 1-2 mjy beam" , for a beam with angular resolution 
of 19'.'2. 126 SMGs were detected at > 3.7a significance (equiv- 
alent to a false-detection rate of ~ 4%,|WeiB et al. 2009) and ro- 
bust or tentative radio, 24 \im or IRAC mid-infrared counterparts 
are id entified to 93 (71 robust and 22 tentative) SMGs (Big gs et"afl 
2010). Here we determine photometric redshifts for the 91 (68 ro- 
bust and 23 tentative) of these SMGs with detectable optical and 
near-infrared counterparts in new and archival multiband photome- 
try of the ECDFS (described in $2^. LESS is an ideal survey for this 
purpose because of its panoramic, deep and uniform submillime- 
tre coverage and extensive auxiliary data, including spectroscopy 
of sufficient SMG counterparts to adequately test our photometric 
redshifts. In addition, the large size of the survey allows us to sta- 
tistically measure the redshift distribution of the SMGs that we are 
unable to locate directly, in order to test if their redshift distribution 
differs significantly from the identified population. 

The plan of the paper is as follows: in ^2] we derive multi- 
band photometry from new and archival observations; while in f]3] 
we describe our the photometric redshift estimates and tests of 
their reliability. The photometric redshifts, SED fits, absolute H- 
band magnitudes, infrared luminosities, dust temperatures and star- 
formation rates of SMGs are presented and discussed in ij4]and we 
present our conclusions in fj5] Th roughout th i s pape r we use de- 
boosted submillimetre fluxes from lWeiB et all f2009b . 12000 coor- 
dinates and ACDM cosmology with Qm ~ 0.3, Q,a = 0.7 and 
Ho — 70kms -1 Mpc~ 1 . All photometry is on the AB magnitude 
system, in which 23.9 iiiab = lp-Jy, unless otherwise stated. 



2 OBSERVATIONS AND DATA REDUCTION 

In this paper we consider the optical and infrar ed counterparts to 
126 SMGs in the ECDFS detected at $ 3.7cr dWeifi et al.ll2009b 
and id entified by VLA ra dio, MIPS dRieke et al.ll2004b 24 urn and 
IRAC dFazio et alj|2004l) emissio n dBiggs et al.ll2010b . Following 
convention and Biggs et alJ J2010t) we consider robust counterparts 
as those with a corrected Poissonian prob ability of being unas so- 
ciated with the submillimetre source (p; iDownes etail 1 19861) of 
p ^ 0.05 in one or more of the radio, 24u.m, or IRAC datasets, 
or p — 0.05 - 0. 10 in two or more; tentative counterparts are those 
with p — 0.05 - 0. 10 in only one of the three bands. 

Six of the SMGs have multiple robust counterparts; of these 
four SMGs (LESS 2, LESS 27, LESS 49 and LESS 74) have two 
counterparts with photometric redshifts ( Sj3. It consistent with them 
being at the same distance and possibly physically associated. Two 
SMGs (LESS 10 and LESS 49) each have two robust counterparts 
with photometric redshifts and SEDs that suggest they are not phys- 
ically associated. In these cases, from the information currently 
available, it is not possible to determine which of the two coun- 
terparts is the source of the submillimetre flux, or whether the 
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Table 1. Summary of photometry employed in this paper. 



Filter 



^effective Detection limit Reference 
(|^m) (3a; mag) 



MUSYC WFI U 0.35 26.9 

MUSYC WFI 1/38 0.37 25.4 

VIMOS U 0.38 28.4^ 

MUSYC WFI B 0.46 26.8~ 

MUSYC WFI V 0.54 26.7 

MUSYC WFI R 0.66 25.8 

MUSYC WFI / 0.87 24.9 

MUSYC Mosaic II z 0.91 24.5 

MUSYC ISPIJ 1.25 23.6 

HAWK-IJ 1.26 25.7 

MUSYC Sofl H 1.66 23.0 

MUSYC ISPI K 2.13 22.7 

HAWK-I if s 2.15 25.3 

SIMPLE IRAC 3.6|J.m 3.58 24.6 

SIMPLE IRAC 4.5|^m 4.53 24.4 

SIMPLE IRAC 5.8^m 5.79 22.8 

SIMPLE IRAC 8.0|t.m 8.05 23.5 



Taylor et al. (2009b) 
Taylor et al. (200%) 
Nonino et al. (2009) 
Taylor et al. (2009b) 
Taylor et al. (2009b) 
Taylor et al. (2009b) 
Taylor et al. (2009b) 
Taylor et al. (2009b) 
Taylor et al. (2009b) 
Zibetti et al. (in prep.) 
Taylor et al. (2009b) 
Taylor et al. (2009b) 
Zibetti et al. (in prep.) 
Damen et al. (2010) 
Damen et al. (2010) 
Damen et al. (2010) 
Damen et al. (2010) 



a The listed depth of the VIMOS U band is that of the central region. The 
typical depth in the shallower outskirts is 28.0 mag 

LABOCA detection is a blend of emission from two galaxies. To 
avoid bias we have included all of the multiple counterparts in our 
analysis, but we note that their small number means that their in- 
clusion does not significantly affect our results. 



2.1 Optical and infrared photometry 

SMGs typically have f aint optical and near-infrared counterparts 
(e.g. llvisonet alJ[2002h so we require deep photometry for ac- 
curate photometric redshift estimates. The ECDFS was chosen 
for this survey because it is an exceptionally well-studied field, 
and as such we are able to utilise data from extensive archival 
imaging and spectroscopic surveys. For completeness and unifor- 
mity we only consider surveys that cover a large fraction of the 
ECDFS rather than the smaller and deeper central CDFS region. 
Therefore, we utilise the MUlt iwavelength Survey by Yale-Chile 
(MUSYC; Gawiser et ai]l2006l) near-infrared survey for U to A- 
band imaging jTavlor et al.l2009bl) . and the Spitzer IRAC/MUS YC 
Publi c Legacy in ECDF S (SIMPLE) imaging for Spitzer IRAC 
data dPamen et alJl2010T) . We also include t/-band data from the 



deep GOODS/VIMOS imaging survey of the CDFS jNonino et all 
l2009h ; although this covers only ~ 60% of LESS SMGs it is valu- 
able for galaxies that are undetected at short wavelengths in the 
shallower MUSYC survey. 

In addition, we have carried out deep near -infrared observa- 
tions in the J and A s bands with HAWK -I JPirard et all 120041 : 
ICasali et alj|200l : iKissler-Patig et alj|2008h at the ESO-VLT (ID: 
082.A-0890, P.I. N. Padilla). The ECDFS was covered with a mo- 
saic of 16 pointings in each band, with a total exposure time of 0.75 
and 1 . 1 hours per pointing, in the J and A s bands respectively. The 
median seeing is 0.7" in J and 0.5" in K 3 . Data reduction has been 
performed using an upgraded version of the official ESO pipeline 
for HAWK-I, customized calibration has been obtained from obser- 
vations of photometric standard stars. More details and catalogues 
will be published in Zibetti et al. (in preparation). 

For accurate photometric redshifts we require consistent pho- 
tometry in apertures which sample the same emitting area in each 



of the 17 filters. For consistency between surveys and to ensure 
that all detected SMG counterparts are included in this study we 
extract photometry from the available survey imaging rather than 
relying on the catalogued sources. SMGs are typically brighter at 
mid-infrared than optical wavelengths due to their high redshifts 
and extreme dust obscu ration. Therefore, we use SEXTRACTOR 
( iBertin & Ar nouts 1996]) to create a source list from a combined 
image of the four IRAC channels, which is weighted such that a 
given magnitude receives equal contributions from all of the in- 
put images. Real sources are required to have at least 4 contiguous 
0.6" x 0.6" pixels with fluxes at least 1.5 times the background 
noise. In addition, we visually check the area within 15" of each 
LABOCA source to ensure that no potential SMG counterparts are 
missed. 

We next use APPHOT in IRAF to measure the fluxes in 3.8" 
diameter apertures for each of the four IRAC bands. We then cut 
the catalogues to 3u based on the background noise, and fi- 
nally apply apertur e corrections as derived by the SWIRE team 
dSurace et alJl2005T) to obtain total source magnitudes. The reso- 
lution in the U- to A'-band imaging is better than IRAC (FWHM 
^ 1.5" compared to ~ 2" for IRAC) and so we convolve each 
U- to A'-band image to match the 1.5" seeing of the worst band. 
We next use APPHOT to measure photometry in 3" diameter aper- 
tures at the positions of the IRAC-selected sources. In all cases, we 
only allow APPHOT to re-centroid the aperture if centroiding does 
not cause the extraction region to be moved to a nearby source, 
as flagged by IRAF's CIER parameter when the centroid shift is 
> 0.5". We have not performed any deblending of the photometry 
but examination of the images suggests fewer than ~ 10% of the 
SMG counterparts are affected. We note here that the photometric 
extraction process is not restricted to SMGs and yields photometry 
(which allows us to calculate consistent photometric redshifts) for 
IRAC-selected sources throughout the ECDFS. 

Finally, to ensure equivalent photometry between the IRAC 
and optical-to-near-infrared filters we create simulated IRAC im- 
ages of point sources. Using these images we calculate that the 
correction between the measured IRAC total magnitudes and the 
photometry extracted from 3" diameter apertures on 1.5" seeing 
images is — 0.014±0.017 magnitudes, and as such we do not apply 
any systematic corrections to the IRAC magnitudes at this stage. In 
i|3,l| we calibrate the photometry prior to photometric redshift cal- 
culation in a process which corrects for small residual offsets. A 
summary of our photometry is presented in TableQ] 

The median number of photometric filters per SMG counter- 
part is 15 and we require detections in at least three photometric 
filters in order to calculate photometric redshifts. Our final sample 
therefore contains 74 optical counterparts to the 68 robustly iden- 
tified SMGs with sufficient detectable optical-to-infrared emission. 
In ij3.3l we show that the exclusion of the additional 23, tentatively 
identified, counterparts does not bias our results. 



2.2 Spectroscopy 

We employ spectroscopy of the ECDFS to calibrate our photometry 
with the SED templates ( i]3.1t and test our photometric redshifts 
( ij3.2t . We have examined th e spectroscopic redsh i ft catalogues 
from many archival surv e ys dCristiani et a l. 2000; Croo m et al. 
| 200l | ; iBunker et all 120031 ; [Dickinson etal] | 2004| ; iLe Fevre et aT 
I2OO4I: 



2004; 



Stan wav et al.l |2004 IStrolger et all |2004|; ISzokolv et al 



van der Wei et al.l |2004|; Izheng et all 12004: iDaddi et al 



2005; Dohertv et al.l2005l;lMignoli et alj200alGrazian et alJ200( 



Ravikumar et al.l 12001 iKriek et alj 120081: IVanzella et alj |20Q: 
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IPopesso et ai1l2009l ; iTreister et alj|2009l ; iBalestra et alj|20ict Ko- 
posov et al. in prep.) and also our own on-going spectroscopic sur- 
vey of LESS sources with the VLT (PID: 183.A-0666, P.I. I. Smail), 
which will be published in full in Danielson et al. (in prep.). 



3 ANALYSIS 



3.1 Photometric redshift calculation 



We use HYPER^I (Bo lzonella et al. I l2000l) to calculat e the photo - 
metric redshifts of counterparts to LESS SMGs dBiggs et al.l2010h . 
HYPERZ compares a model SED to observed magnitudes and com- 
putes x 2 f° r eacn combination of spectral type, age, reddening 
and redshift and thus statistically determines the most likely red- 
shift of the galaxy. We use the elliptical (E), Sb, single burst 
(Burst) and constant star -formation (Im) spectral templates from 
iBruzual & Chariot d 19931) which are p rovided with HYPERZ, and 
allow reddening 1 Calzetti e t al. 2000) of Ay — 0-5 in steps 
of 0.2. This combin ation of templates and Ay was shown by 
Wardl ow et alj d2010h to be sufficient for calculating photometric 
redshifts of SMGs. Redshifts between and 7 are considered and 
galaxy ages are required to be less than the age of the Universe at 
the appropriate redshift. In i]3.5| we show that the HYPERZ-derived 
galaxy ages cannot be reliably determined, but we note here that 
the requirement for SMGs to be younger than the Universe does not 
significantly affect the derived redshifts. Galaxies are assigned zero 
flux in any filter in which they are not detected, with an error equal 
to the lcr detection limit of that filter. To ensure that galaxies at 
z ~ 2-3 do not have their redshifts systematically underestimated 
we have modified the handling of the Lyman-a forest in HYPERZ, 
such that intragalactic absorption in the models is increased and 
three different levels of absorption are considered in the fitting pro- 
cess. The reliability of the calculated redshifts and the validity of 
these settings is tested in full in €|3 .21 

We test for small systematic discrepancies between the pho- 
tometry and model SEDs prior to using HYPERZ to calculate pho- 
tometric redshifts of SMGs. This is done by running HYPERZ on 
1796 galaxies and AGN with spectroscopic redshifts in the ECDFS 
and requiring a fit at the observed redshift. We then compare the 
model and measured magnitudes for each galaxy, and iteratively 
adjust the zeropoints of the filters with the largest systematic off- 
sets. This yields significant offsets for the following filters: VIMOS 
U (0.083 mag), MUSYC U (-0.091 mag), C/38 (-0.074 mag), R 
(0.049 mag), / (0.048 mag), z (0.095 mag), HAWK-I J (0.043 
mag), IRAC 3.6 (0.043 mag) and IRAC 8.0 p.m (0.110 mag). The 
typical uncertainties in these corrections are ±0.02 and the remain- 
ing eight filters have no significant corrections. 

The calibrated photometry of the robust LESS SMG counter- 
parts is listed in Tableland in Table|2]we provide the coordinates, 
photometric redshifts, absolute rest-frame //-band magnitudes, far- 
infrared luminosities and characteristic dust temperatures of the 
SMGs ( fl4~6l . We also provide the reduced \ 2 of the best fit SED at 
the derived photometric redshift and the number of filters in which 
the SMG was detected and undetected (but observed). We caution 
that the reduced x 2 for galaxies with only a few photometric de- 
tections is typically low (< 0.5) but the error on the photometric 
redshift is typically large, since there are only weak limits on the 
SED from the photometry. Therefore, the values of the reduced x 2 
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Figure 1. Spectroscopic redshift against Az/(l+z) for robust counterparts 
to LESS SMGs. We distinguish between high and low quality spectroscopic 
redshifts as determined by the flags provided in most archival catalogues, 
and highlight the two likely quasars (LESS 66 and LESS 96; see Appendix 
|7\}. The median (mean) Az/(1 + z) for the all the SMGs is 0.023 ± 0.021 
(-0.013 ± 0.178). The inset plot shows the histogram of Az/(1 + z) 
for 1796 galaxies and AGN in the ECDFS with spectroscopic redshifts. 
The distribution is centered on 0.016 ± 0.002 and has a la dispersion 
of 0.05. We conclude that our photometric redshifts are a good proxy for 
spectroscopic redshifts for both samples. 



should be considered in conjunction with the number of photomet- 
ric detections when considering the reliability of the photometric 
redshifts. 

The median reduced x 2 of the SMG counterparts is 2.3 (2.1 
if only the galaxies with reduced x 2 ^ 10 are considered). This 
suggests that our photometric errors are slightly overestimated and 
lead to apparently overly-precise photometric redshift limits. In- 
deed, we find that the HYPERZ 99% confidence intervals more 
reliably represent the la errors, yielding ~ 68% of SMGs with 
photometric redshifts consistent with the spectroscopic redshifts. 
Therefore, throughout this paper we use the HYPERZ 99% confi- 
dence intervals on the photometric redshift estimates to represent 
the 1-er uncertainty. Of the 74 SMG counterparts examined there 
are eight with poor fits of the SED to the photometry (indicated 
with reduced x 2 > 10). Of these, one (LESS 39) is blended in the 
optical imaging and two (LESS 66 and LESS 8 1 ) lie in stellar halos. 
LESS 66 is also likely to be a quasar, as is LESS 96, and another 
four SMGs with reduced X 2 > 10 (LESS 19, LESS 57, LESS 75 
and LESS 111) have excess 8 urn flux compared to the best-fit SED, 
which is indicative of an AGN component (see ^3.4l for a full dis- 
cussion). Since we did not include any quasar or AGN templates 
in the fitting procedure it is unsurprising that these sources are not 
well represented by the employed SEDs. We note here that, as we 
show in ^3.41 the exclusion of AGN templates does not bias our 
photometric redshift estimates. 



3.2 Reliability of photometric redshifts 

To test the reliability of our photometric redshifts (z p hot) we first 
compare them to the spectroscopic redshifts (z S p CC ) for 1796 galax- 
ies and AGN in the ECDFS and calculate Az — z spac — z p hot for 
each source. The histogram of Az/(l + z) for these 1796 sources is 
shown as an inset in FigJTJ the sample is centered on Az/(l + z) = 
0.016 ± 0.002 and has a la dispersion of 0.05. We define outliers 
as sources with | Az|/(1 + z) > 0.3; the outlier fraction for these 
1796 field galaxies and AGN is 0.15. We also calculate the outlier 
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Table 2. The catalogue of 74 robust counterparts to LESS SMGs, their photometric redshift estimates, reduced \ 2 of the best-fit SED and the number of 
photometric filters in which the galaxy is observed. We also present the absolute rest-frame H-band magnitudes, the derived far-infrared luminosities and 
characteristic dust temperatures of the SMGs. 
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-24. 


.04 


LESSJ033248 


.1- 


■275414 


LESS 12 


03 h 32 m 47!96 


-27°54'16'.'l 


3 


„,^4-1 


112 
11 


0.1 


6 [11] 


-24. 


.06 


LESSJ033152 


.6 


■280320 


LESS 14 


03 h 31 m 52?47 


-28°03'18'.'6 


3 


' 56 ±0 


<!_> 
.")(> 


0.8 


7 [9] 
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-0.4 
^+4.4 

-5.4 

10.0II 8 

11.91*:? 

1+7.9 
-3.7 

l 2 -3le 2 2 2 

0+2.2 
5 -1.8 
1+1.6 
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34.5^ 
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58.4 
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4.3 
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43.6 
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31.8 
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23.8 
48.3 
41.0 
42.8 
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49.3 
29.2 
29.6 
43.1 
25.4 
54.4 



+3.6 
3.5 
+3.6 
2.5 
+2.7 
3.8 
3.1 
4.2 
+2.4 
2.4 
+5.4 
0.1 
+3.4 
3.4 
+3.2 
3.4 
+5.2 
4.6 
+9.9 
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+4.3 
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.3.4 
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+3.4 
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+5.0 
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Table 2 - continued 



SMG a 



Short name 



RA" 



Dec 6 



^phot 



2 d 
A rcd 



Filters 6 



(mag) 



Lfir 

12 



(10 



a 



T D h 
(K) 



ID type 1 



LESSJ033154.2- 
LESSJ033251.1- 
LESSJ033 155.2- 
LESSJ033313.0- 
LESSJ033130.2- 
LESSJ033151.5- 
LESSJ033335.6- 
LESSJ033325.4- 
LESSJ033 140.1- 
LESSJ033316.4- 
LESSJ033122.6- 
LESSJ033325.6- 
LESSJ033249.3- 
LESSJ033 150.8- 
LESSJ033 128.0- 
LESSJ033121.8- 
LESSJ033328.5- 
LESSJ033139.6- 
LESSJ033209.8- 



275109 
273143 
275345 
275556 
275726 
274552 
274020 
273400 
275631 
275033 
275417 
273423 
273112 
274438 
273925 
274936 
275655 
274120 
274102 



LESS 84 
LESS 87 
LESS 88 
LESS 96 
LESS 98 
LESS 101 
LESS 102 
LESS 103 
LESS 106 
LESS 108 
LESS 110 
LESS 111 
LESS 112 
LESS 114 
LESS 117 
LESS 118 
LESS 120 
LESS 122 
LESS 126 



03 h 31 m 54!49 
03 h 32 m 50f83 
03 h 31 m 54?81 
03 h 33 m 12!62 
03 h 31 m 29?89 
03 h 31 m 51?53 
03 h 33 m 35!56 
03 h 33 m 25?37 
03 h 31 m 40fl7 - 



03 n 33' 
03 h 3T 



M6?51 
J 22!63 



03 h 33 m 25!21 
03 h 32 m 48?85 
03 h 31 m 51?08 
03 h 31 m 27!62 
03 h 31 m 21!91 
03 h 33 m 28?55 
03 h 31 m 39?52 
03 h 32 m 09?60 



27°51'05'/3 
27°31'4h"2 
27°53'40'.'9 
27°55'51"6 
27°57'22"4 
27°45'53('l 
27°40'23y2 
27°33'58'.'5 
27°56'22'.'4 
27°50'39 / /3 
27°54'17"0 
27° 34' 25^9 
27°31'12?8 
27°WiT.'0 
27°39'27"3 
27°49'34"0 
27°56'54"1 
27°41'19^'4 
27°41'06?9 



2.29 
3.20 
2.35 
2.71 
1.55 
2.39 
1.68 
1.84 
1.96 
0.20 
2.35 
2.61 
1.81 
1.57 
1.73 
2.17 
1.43 
2.08 
2.02 
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-0. 
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.+0. 
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,+0. 


29 


'-0. 


34 


■+4. 


S3 




.49 


■+<)■ 


.30 


'-0. 


21 


,+0. 


us 


'-0. 


OS 


,+o. 


17 


'-0. 


13 



3.6 
0.1 
1.1 
22.0 
1.0 
2.5 
1.1 
0.3 
2.1 
6.3 
0.0 
14.4 
0.7 
1.6 
3.3 
1.7 
2.2 
5.2 
2.4 



14 [3] 
5 [0] 

16 [1] 

17 [0] 
10 [4] 

10 [7] 

11 [2] 
5 [7] 

11 [5] 

15 [0] 
4[0] 
13 [0] 
5 [0] 
17 [0] 
9 [4] 
5 [1] 
13 [3] 
17 [0] 

12 [4] 



-24.14 
-24.84 
-24.37 
-26.30 
-24.40 
-23.51 
-24.34 
-23.44 
-25.00 
-22.75 
-23.22 
-24.49 
-24.02 
-24.61 
-24.23 
-22.21 
-23.41 
-25.14 
-23.84 



+ 1.8 

1.2 
+ 5.6 
19.3 
2.3 
-1.9 
+ 1.9 
-2.5 
2.S 
2.1 
2.9 
2.3 



4.5 
37.0 
11.7 
16.0 

7.8 

3.8 

< 1.3 

< 1.6 

6 3+ 34 
°- d -3.4 
n 9 +0.08 
u - -0.08 

< 2.8 
9.8 
2.3 
5.3 
5.7 
2.8. 
2.1 
22.4 
2.3 



60.1 
44.0 



+3.4 
2.3 
+ 2.5 
1.2 
1.1 
0.9 
+3.4 
2.6 
+47.5 
2.7 
-1.8 
-0.9 
+2.8 
2.7 
+ 1.3 
0.9 



34.5±|;S 

+5.9 

-12.9 
+4.0 

49.7:!| 
39.9±t;i 
33-8±l:| 

< 24.9 

< 26.3 
38.5^4 

24.5±^3 

< 31.4 

5.0 
4.S 



44.1 
28.5 
36 
37.7 
31.8 
29.2 
55.2 
30.6 



4.9 
+3.7 
-3.7 
+5.5 
6.0 
+48.7 
-15.6 
+5.0 
4.3 
+6.2 
6.2 
+4.2 
4.1 



I 

RM 
R 

RM 

RM 

R 

M 

M 

RI 

RM 

MI 

RM 

RI 

RM 

R 

R 

RM 
RM 
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"The SMG names correspond to those in lWeiB et al] <2009h and lBiggs etafl feplOh . 
^Coordinates are the J2000 position of the optical/near-infrared counterpart. 

c Since HYPERZ was restricted to < z < 7 the six galaxies whose upper redshift limits yield a formal maximum redshift of z max = 7 are actually only 
constrained in the lower redshift limit. Therefore, throughout this paper the redshifts of these galaxies are plotted as lower limits. 
d The reduced \' 2 of the best-fit SED at the derived photometric redshift. 

e The number of photometric filters in which each SMG counterpart was detected [and the number of filters in which the SMG was observed but not 

detected, providing a limiting flux] . 

* M h is the absolute magnitude in the rest-frame //-band. 

9 As discussed in £|4.6l the far-infrared luminosity (Lfir) is derived from the infrared-radio correlation using the radio flux and the photometric redshift of 
each SMG. 

h The characteristic dust temperature (To) is derived as discussed in ij4.6l from radio and sub millimetre fluxes and the photometric redshift of each SMG. 
"ID types R, M and I indicate radio, 24 i^m and IRAC identified counterparts respectively (see iBiggs et aljioTol for details). 

■'As shown in £|4. 1I LESS 20 appears to contain a radio-loud AGN. Therefore, the Lpm and To presented here are likely significantly overestimated due 
to the AGN contribution to the radio flux, as such LESS 20 is excluded from our studies of the luminosity function, star-formation rates and star-formation 
history of SMGs ( jgg) , 



resistant normalised median absolute deviation (NMAD) of Az, 
(tnmad = 1.48xmedian(|A2— median(Aj:)|/(l+2)) = 0.097. 
These statistics show that our photometric redshifts are a good 
proxy for spectroscopic redshift for these sources. However, the 
median redshift, z = 0.84, is lower than that expected for SMGs 
and the targets are typically brighter at optical wavelengths, limit- 
ing the usefulness of these comparisons for the SMGs. Therefore, 
we also test our photometric redshift calculation on the 30 robust 
SMG counterparts with available spectroscopic redshifts. 

Fig. [T] shows spectroscopic redshift against Az/(1 + z) for 
the 30 robust LESS SMG counterparts with spectroscopic redshifts 
from archival surveys of the ECDFS (10 counterparts) and our 
spectroscopic survey of SMGs in the ECDFS (20 counterparts). 
Quality flags are published in many catalogues and where possi- 
ble we distinguish between high- and low-quality spectroscopic 
redshifts. The median Az/(1 + z) for SMGs is 0.023 ± 0.021 
(here are throughout this paper errors on median measurements are 
from bootstrapping) the mean Az/(1 + z) = —0.013 ± 0.178, 
and (tnmad = 0.037, suggesting that for SMGs our photometric 
redshifts are reliable. We caution that the SMGs without reliable 



spectroscopic redshifts are fainter on average than the SMGs with 
spectroscopic redshifts, which could affect the quality of their pho- 
tometric redshifts. Although the median .R-band magnitude of the 
SMGs detected in the MUSYC survey is the same for the counter- 
parts with and without spectroscopic redshifts, all of the 30 SMGs 
with spectroscopic redshifts are detected in the MUSYC R band, 
while only 18 of the 45 SMGs without spectroscopic redshifts are 
detected. 

|PunlopetalJ ( l2009l) have independently calculated photomet- 
ric redshifts for the six LESS SMGs (five with robust counter- 
parts) in GOODS-South th at were also detected by BLAST at 
250 (im dPevlin et al]|2009h . Their photometry uses imaging from 
HST (B435, Veoe, «775, zsso), the VLT (J, H, K) and Spitzer (3.6, 
4.5, 5.8 and 8 |i.m) and they use HYPE RZ with the st ellar popula- 
tion models of Chariot & Bruzua l (e.g. Bruzu al2007t) which have 
a Salpeter IMF. In all cases the (Du nlop et al.ll2009l) photometric 
redshifts agree with those presented in Table [2] providing further 
confidence that our photometric redshifts are reasonable for SMGs. 

To assess the level of systematic uncertainties for the derived 
photometric redshifts due to the adopted methodology, SED tem- 
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plates, and/or photometric data, we also u se the Zurich Ex tragalac- 
tic Bayesian Redshift Analyzer (ZEBRA; lFeldmann et~aT]|2006l) to 
calculate photometric reds hifts. Our adopted procedure is similar 
to that discussed in §3.2 of lLuoetalJ feoid) . Briefly, we use ZE- 
BRA to obtain a maximum-likelihood estimate for the photomet- 
ric redshifts of individual galaxies or AGNs using an initial set of 
265 galaxy, AGN, and galaxy/ AGN hybrid SED templates. These 
SED templates were then expanded to 463 templates during the 
template-training mode of ZEB RA to best repre sent the SEDs of 
the « 2 Ms CDFS X-ray sources dLuo et al.l2 010). including AGN. 
Besides the different SED templates used, this method differs from 
the HYPERZ approach described in t]3.1| in some additional details 
such as how the redshift int ervals and minim um photometric errors 
are determined; see §3.2 of iLuo et"al I d2010h for details. 

The ZEBRA-derived photometric redshifts (zphot, check) were 
compared to those listed in Table 1 (.Zphot); the difference was 
measured by <5z ph ot = (2 P hot,check - z P hot)/(l + ■^phot 

). For 

sources with secure spectroscopic redshifts, individual |<52 p hot| val- 
ues range from as 0.01 to 0.10, indicating that both methods are 
able to deliver photometric redshifts to a similar accuracy. For the 
full sample, the mean (median) value of <5z p hot is —0.006 (0.011), 
with an rms scatter of 0.028, suggesting that the photometric red- 
shifts in Table 1 are fairly robust. After accounting for the effective 
la errors of the photometric redshifts, only three (sources LESS 7, 
LESS 37 and LESS 1 1 1) of the 74 sources have inconsistent z p hot 
and %,hot, check- As some sources have photometry data in addition 
to those presented in Table 2, we also tested the effect of includ- 
ing more data points. The photometric redshifts differ by a mean 
value of |<5z p hot| of 0.024 with an rms error of 0.030, after includ- 
ing th e WFI 7?-band data dGiacconi et al.l |2002| ; iGiavalisco et al.l 
120041) for 25 sources and the GALEX near-UV and far-UV data 
1 Morr issev et al.l F2007) for 3 sources. Given the small difference 
caused by the additional data, we consider the consistent aperture 
photometry in Table 2 suitable for the purpose of deriving reliable 
photometric redshifts. 



3.3 Sample subsets 

As discussed in SjT] Biggs et al.l J2() 1 (■ ) identified robust counter- 
parts to 71 SMGs (of which 68 have detectable optical counter- 
parts) and tentative counterparts to 26 LESS SMGs, respectively. 
In Fig. [2] we compare the redshift distributions of the robust and 
tentative counterparts to determine whether our results may be bi- 
ased by the exclusion of tentative counterparts in our main anal- 
ysis. The median redshifts of robust and tentative counterparts are 
2.2±0.1 and 2.0±0.2 and the interquartile ranges are z = 1.8-2.7 
and z = 1.2-2.8, respectively. We also use a Kolmogorov-Smirnov 
test (KS-test), which calculates the probability that two samples are 
drawn from the same parent population (Pks), to compare the red- 
shift distributions of robust and tentative counterparts statistically. 
We find Pks = 0.09 and conclude that it is likely that there is some 
contamination from physically unassociated foreground (z < 1) 
galaxies in the tentative identifications. Therefore, throughout the 
remainder of this paper we restrict our analysis to robust counter- 
parts. 

We note that our identified sample may contain a small 
number of potential gravitational lenses. These are typically low- 
redshift counterparts where the radio or mid-infrared emission is 
offset from the optical source. These are discussed individually in 
Appendix [A] and we have confirmed that their inclusion does not 
affect our results. 

In Fig. [2] we also compare the redshift distributions of coun- 



Field galaxies (rescaled) 
Robust IDs 
Tentative IDs 




Figure 2. In the top panel we present a comparison between robustly 
and tentatively identified SMG counterparts and the field population of 
the ECDFS; tentative counterparts and field galaxies are offset slightly 
in redshift for clarity. The robust counterparts have a median redshift of 
z = 2.2 ± 0.1, compared to z = 2.0±0.2 for the tentative counterparts. 
Tentative counterparts have a larger interquartile range - 1.2-2.8 compared 
to 1.8-2.7 for the robust SMG counterparts. We interpret these distributions 
as evidence that tentative counterparts are mainly drawn from the same par- 
ent population as the robust counterparts, but with the addition of some 
contamination, particularly at low redshifts. In the lower panel we compare 
the redshift distributions of robust counterparts (with p ^ 0.05) in radio, 
24 |xm and IRAC data; for clarity radio and 24 u.m counterparts are plotted 
offset slightly in redshift. The radio, 24 |j.m and IRAC samples have median 
redshifts of 2.3 ± 0.1, 2.1 ± 0.2 and 2.3 ± 0.2 and interquartile ranges of 
1.8-2.8, 1.6-2.5 and 1.9-2.7 respectively. Therefore, we find no significant 
differences in the redshift distributions of the three identification methods. 



terparts with p ^ 0.05 in the radio, 24 \im and IRAC data. The 
median redshifts are z = 2.3 ± 0.1, 2.1 ± 0.2 and 2.3 ± 0.2, and 
the interquartile ranges are z = 1.8-2.8, 1.6-2.5 and 1.9-2.7 for 
the radio, 24|i.m and IRAC samples respectively. A comparison of 
the three redshift samples shows that they are statistically indis- 
tinguishable. We conclude that the three counterpart identification 
methods select galaxies with similar redshift distributions, and are 
not significantly biased with respect to each other. 



3.4 The effect of AGN on photometric redshifts 

As discussed in &13. 1 l our photometric redshift calculations are based 
on fitting stellar templates to the SMG photometry. However, stud- 
ies have shown that the 8 vim flux in SMGs with a luminous AGN 
component can be dominated by the AGN an d therefore fitting stel- 
lar templates may yield m isleading results l lHainline et at] [2009, 
l2010l ; ICoppin et al.l(2010bl) . 

We employ two methods to identify potential AGN in the 
LESS SMGs. Firstly, we cross-correlate the LESS SM G counter- 
parts with the Cha ndra X-ray catalogu es of the CDFS dLuo et all 
2008) and ECDFS dLehmer et alj F2005) with a matching radius of 
1". This yields 12 X-ray luminous SMG counterparts (16% of the 
robust LESS counterparts). Secondly, we identify nine SMG coun- 
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terparts (12% of the robust LESS counterparts) with a large excess 
of 8 \ua flux compared to the best-fit SED template, which poten- 
tially indicates obscured power-law emission from an AGN. In total 
this yields 16 SMG counterparts (22% of the robust LESS counter- 
parts) which may contain AGN (five are both X-ray detected and 
have an 8 [im excess). 

To determine whether AGN contamination in the S-[im filter 
reduces the accuracy of our photometric redshift estimates we re-fit 
the photometry of the counterparts which display an 8-iim excess 
whilst excluding the 8 iim photometry. Spectroscopic redshifts are 
available for eight of the affected galaxies and we find that the av- 
erage \Az\ of these eight galaxies is not significantly decreased, 
while the median reduced \ 2 drops by 1.2, when the 8 vim photom- 
etry is excluded from the fitting. This result does not change if we 
also exclude the 8 vim photometry of the SMGs which are X-ray 
detected. 

It is also possible that our sample of SMG counterparts con- 
tains AGN which enhance the 8 vim flux but do not cause a de- 
tectable excess. Therefore, we also exclude the 8 vim photometry of 
all the SMG counterparts during the fitting procedure. Once again 
the average \Az \ for the spectroscopic sample does not change sig- 
nificantly (median Az/(1 + z) = 0.032 ± 0.021, compared to 
median Az/(1 + z) = 0.023 ± 0.021 originally). These results 
indicate that when calculating photometric redshifts the benefit of 
including the longer-wavelength data is greater than the bias which 
is removed by ignoring the 8 vim photometry. Therefore, we include 
the 8 urn photometry in the SED fitting. 



3.5 Reliability of SED parameters 

In addition to calculating photometric redshifts, HYPERZ also re- 
turns the spectral type, age and reddening of the best-fit SED tem- 
plate. To test the sensitivity of the choice of template we refit the 
photometry of the SMGs allowing only the Burst template, and then 
only a constant star-formation rate history (Im). These two tem- 
plates represent the extremes of the star-formation histories and so 
they will let us gauge the sensitivity of the derived parameters to 
the choice of the best-fit template. 

We compare the quality of the Burst and Im fits of each galaxy 
with AXrcd ~ m e difference in the reduced \ 2 of the Burst and Im 
fits. We find that 63% of the SMGs have | AxLiI ^ 1 and as such 
the Burst and Im templates are indistinguishable at the 99% level 
for these SMGs. 61% (17) of the SMGs with |Ax? c d| > 1 between 
the two template fits are best-fit by Bursts and 39% (11) by Im 
templates. If three templates - Burst, Sb and Im - are considered 
the star-formation histories of only 23% (17) of all the SMG coun- 
terparts can be distinguished. Therefore, although it may be pos- 
sible to crudely distinguish the star-formation histories of a frac- 
tion (~ 20^10%) of SMG counterparts with HYPERZ, the star- 
formation histories of most SMGs counterparts cannot be reliably 
established. 

We find that the SMGs that have SED fits with \A X 2 cd \ < 1 
and can be equally well fit by either Burst or Im templates have dif- 
ferent age estimates depending on the template. The age and star- 
formation history of a stellar population affects the light-to-mass 
ratio; thus our inability to distinguish between star-formation his- 
tories leads to uncertainties in the light-to-mass ratios and stellar 
mass estimates. Using the ff-band light-to-mass ratios for Burst 
and Im models from the STARBURST99 stellar population model 
dLeitherer et"ai]|l999h we calculate that the uncertainties in SED 
fitting parameters result in a la dispersion of a factor of 4.6 range 



in the light-to-mass ratios and consequently in the stellar mass es- 
timates (see also Fig. [9}. 

We also compare reddening measurements for those galaxies 
with SED fits with |Ax 2 | ^ 1 and find that on average the differ- 
ence between Ay for the best-fit Im and Burst templates, A Ay, 
is equal to 0.32 ± 0.16. Since average estimates based on either 
template return the same value of Av (to ~ 2a) we conclude 
that average reddening measurements for the SMG population are 
not strongly sensitive to the adopted star-formation history and are 
likely to be statistically meaningful (although we caution against 
trusting values for individual SMGs). When fitting E, Sb, Burst, 
or Im templates and allowing Av ~ 0-5, as in our photometric 
redshift calculations we determine a median Av = 1.5 ± 0.1 and 
96% of SMGs have Av ^ 3. We conclude that in most instances 
restricting to Av ^ 3 is sufficient for calculating photometric red- 
shifts of SMGs. We note that these values of Av are integrated 
across the whole galaxy and that obscuration in the regions respon- 
sible for the majority of the far-infrared/submillimetre emission is 
considerably larger (e.g. IChapman et al]|2004bl : iTakata et alj|2006l ; 
llvison etalfeOlObl) . 



4 RESULTS AND DISCUSSION 

In ij3]we derived reliable photometric redshifts for 74 counterparts 
to 68 robustly identified SMGs. We now use these photometric red- 
shifts and the SED fits to further investigate the properties of SMGs. 

4.1 Photometric redshifts 

In Fig. [3] we show our photometric redshift distribution for the 74 
robust SMG counterparts; it peaks at z = 2.2 ±0.1 and has an 
interquartile range of 1.8-2.7. We compare to the photometric red- 
shift distribution of SMG counter parts in the SCUBA Half-Degree 
Extragalactic Survey (SHADES) ^Clements et alj|2008l ; iDve et all 
2008), median z = 1 .5 ± 0.1, and the spectroscopic sample from 
Chapman et al. (2005), median z = 2.2 ± 0.1, both of which have 
similar submillimetre flux limits as ou r survey. The L ESS SMGs 
have a similar redshift distribution to Chapman et al. | d2005h . al- 
though in LESS the spectroscopic 'redshift desert' at z ~ 1.2-1.8 
is filled and there is a larg er high-redshif t tail. T he redshift distri- 
butions of both LESS and Chap man et alj {2005) SMGs are peaked 
at higher redsh i fts tha n the SHADES SMGs. A KS-test between 
IChapmanetafl {2005) and LESS SMGs yields P KS = 0.44, sug- 
gesting the two samples appear to be drawn from the same parent 
population. However, a KS-test between the LESS and SHADES 
SMGs gives Pks = 1-3 x 10~ 5 indicating that these samples are 
likely drawn from intrinsically different populations. 

We conclude that the global properties of our photometric red- 
shifts are consistent with the largest previous spectroscopic sur- 
vey, albeit with a higher-redshift tail - we find 10 (14%) SMGs 
with z ^ 3 and eight (11%) with z ^ 3.5, o f which LESS 73 is 
spectro scopically confirmed at z = 4.76 (see Coppi rTet alj|2009l . 
l2010al) . It is likely that the larg er number of high-reds hift sources in 
the LESS survey compared to lChapman et alj d2005h . where there 
are 10% at z ^ 3 and just 1% at z ^ 3.5, is due to deeper ra- 
dio da ta (on average) and the inclusion of 24- vim counterparts in 
LESS (Big gs et al JI2010T) and most critically the use of photomet- 
ric redshifts covering the UV to mid-infrared which are less reliant 
on the detection of spectral features in the optical. Conversely, the 
SHADES SMGs appear to typically lie at lower redshifts. We stress 
that compared to the SHADES analyses we have used about twice 
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Figure 3. [left] The photometric redshift distri bution of the LESS SMGs. We compare this to the photometric reds hift distribution of SHADES SMG coun- 
terp arts IClements e t al. 2008; Dve et al. 2008]), and the spectroscopic redshift distribution of SMGs from lChapman et alj 120051) : for clarity the SHADES 
and lChapmarTetal J |2005|) samples are offset slightly in redshift. The median redshift of identified SMGs in LESS is z = 2.2 ± 0.1, which is the same as 
that from the lChapmanTtal] l2005h spect roscopic survey ; SHAD ES has a lower median redshift of z = 1.5 ± 0.1. There is a slightly larger high-redshift 
tail in the LESS SMG popul ation than the[chapman et al. 12005) SMG population. Additionally, the so-c alled 'redshift desert' at z ~ 1.5, which is evident 
in the lChapman et all j2005h study does not affect our photometric redshifts and as such, in contrast withlChapman et al l |2005D, the increase in the number 
of galaxies from z ~ 1 to z ~ 2 is smooth. Statistical comparisons show that the Chapman et al. (2005) and LESS SMGs are most likely drawn from pop- 
ulations with similar redshift distributions, but that the SHADES SMGs are biased to low redshifts either from systematic errors in the photometric redshift 
calculations, sample selection, or cosmic variance. The shaded region represents the area that would be added to the histogram were the redshifts of the 57 
statistically identified or completely unidentified SMGs known and is designed to give an impression of the potential contribution of the unidentified SMGs 
to this figure. The lower region corresponds to the unidentified SMGs that we statistically identify in £14.21 and which have redshifts similar to the identified 
SMGs. The upper shaded region represents the SMGs which remain unaccounted for after the statistical analysis and likely have z > 3. [right] The same 
distribution but now including the statistically identified SMG population from £14.21 in addition to the identified sample. These are distributed uniformly 
within the relevant Az = 1 ranges. The shaded area now represents the remaining unidentified SMGs, which are likely to lie at z > 3. We conclude that the 
median redshift of the Ss70fj,m > 4mJy SMG population is likely to be z = 2.5 ± 0.6. 



as many photometric bands, tested against a larger spectroscopic 
sample of SMGs, and obtained qualitatively better fits to the SEDs. 
We suggest that either there is a systematic error in the original 
SHADES photometric redshifts or their counterpart identifications, 
or that cosmic variance is the cause of the different redshift distribu- 
tions. However, we note that a re-analysis of the optical-to-infrared 
photometry of the SHADES SMGs yields a median photometric 
red shift of z = 2.05 (Scha el et al. in prep.) - more similar to LESS 
and Chapma n et al.1 d2005h than the original SHADES analyses. 

Studies have suggested that the brightest SMGs may have 
higher redshifts than those with lower submillimetre fluxes (e.g. 
Ilvisonetal .1120021 ; IPope etalj|20ol iBiggs et al.|[20ld) . In Fig. 0] 
we plot the photometric redshift against 870 [im flux (SWo^™) f° r 
robust LESS SMG counterparts. We split the galaxies into those 
brighter and fainter than the median deboosted submillimetre flux 
of the sample, SWo/™ = 5.6 mjy (we use deboosted 870-/im 
fluxes through out). SMGs with Sg7o M m ^ 5.6 mjy have a median 
redshift of z = 2.1 ± 0.2 and SMGs with SWo^m > 5.6 mJy have 
a median redshift of z = 2.3 ± 0.2, where the errors are bootstrap 
uncertainties on the medians. Spearman's rank correlation coeffi- 
cient between SWopm and 2 p h t is 0.20, which corresponds to a 
probability of zero correlation of 0.08 and indicates that there is no 
significant correlation between submillimetre flux and redshift for 
SMGs in our sample. We have verified that the result is not depen- 
dent on the choice of the flux limit between the two bins. Addition- 
ally, if all the unidentified SMGs lie at z — 5 or z — 1 (in 34.21 
both of these scenarios are shown to be unlikely) we still find no 
statistically significant difference between the redshifts of SMGs 
in the two flux bins. The sample of SMGs with optical-infrared 
photometric redshifts in this work is larger than previous studies of 
this phenomenon and our analysis finds no significant correlation 




3 4 5 6 7 8 9 1011 1213 

s 87 o ( m Jy) 



Figure 4. Photometric redshift versus submillimetre flux for LESS SMGs; 
the median 5g70fim error bar is shown in the bottom right. The median 
Ss70iJ.m and redshift, with lcr error bars, are presented for SMGs with 
Ss70iJ,m ^ 5.6 mJy and Ssr0u.m > 5.6 mJy. Previous studies have sug- 
gested that the brightest SMGs may lie at the highest redshifts. This work 
contains optical-infrared photometric redshifts for a larger sample of SMGs 
than previous studies of the phenomenon and finds no evidence for a trend. 
For comparison we also highlight the 1-cr d istribution o f SMGs in fl ux bins 
of 1 mjy in the ACDM GA LFORM model iBaugh et aT]|2005l ; lEacev et alj 
l2008tlSwinbank et alj2008l) . 

between Ssio^m and redshift for robustly identified sources, also 
implying th at Ssyoum is no t a go od proxy for redshift. This result 
agrees with Knu dsen et al.1 fcOlCh . who find no difference in the 
redshift distributions of faint lensed SMGs (Sasoum < 2 mjy) an d 
the brighter (Saso^m > 3 mjy) SMGs from Cha pman et~al1 ( l2005l) . 
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Figure 5. The top panel shows the redshift histogram of sources within the submillimetre positional error circles of SMGs without robust radio, 24 |xm or 
IRAC counterparts compared to the same number of random positions in the field. For comparison we also plot the redshift histogram for sources in the 
submillimetre positional error circles of SMGs with robust counterparts, also scaled to the same number of error circles. In the bottom panel we show the 
difference between the redshift histogram of galaxies around unidentified SMGs and the field population (from random positions). We also plot both the 
difference in redshift of g alaxies aroun d identified SMGs and the field population, and the redshift distribution of radio-undetected SMGs in the ACDM 
GALFORM model JSwinbank et alfcOQa). In order to highlight potential differences in the redshift distributions of the populations the latter two datasets are 
scaled to match the value of AiV ga i of unidentified SMGs in the z = 2-3 bin. By using AN ga \ and assuming a uniform distribution of galaxies within the 
bins we calculate that the average redshift of unidentified SMGs is z = 2.5 ± 0.3. 



We note that SMG s in the semi-analyt i c ACDM GALFORM m odel 
l lBaugh et al]|2005l : lLacev et alJkOOa ISwinbank et al.ll2008h also 
show no correlation between Ssro^m and redshift (although the er- 
ror range decreases at high fluxes where there are few galaxies in 
the model). 

4.2 Redshift distribution of unidentified SMGs 

To date, redshift surveys of SMGs have focused on the ~ 60-80% 
of the population with counterparts identified from radio and 24 |j.m 
imaging, and a few located using high-resolution (sub-)millimetre 
interferometry. The requirement for radio or infrared counterparts 
to SMGs can bias the redshift or the dust te mperature distributions 
of identified SMGs ( Chapman e t alj|200"3) and it is currently un- 
known if the identified population is representative of the ~ 20- 
40% of SMGs without identified counterparts. In particular, it is 
unclear whether they have the same redshift distribution. In order 
to investigate the redshift distribution of the unidentified SMGs we 
utilise our extensive 17-band photometric redshifts in the ECDFS 



to investigate the photometric redshifts of sources around SMGs 
without robustly identified counterparts to the field population. 

In Fig. [5] we show the redshift histogram of galaxies within 
the error circles of the 55 unidentified SMG s in the ECDF S, where 
the region considered is that used by Big gs et al] d2010t) to iden- 
tify SMG counterparts (both the completely unidentified and those 
with only tentative identifications). For comparison, we show the 
redshift histogram of all the galaxies in the submillimetre error cir- 
cles of SMGs with robustly identified counterparts, scaled such that 
the number of error circles examined is the same as the unidentified 
SMG sample. We also consider the photometric redshifts of galax- 
ies in the same area around random positions in the field, which 
are required to be > 15" from any LESS SMGs. We consider 50 
Monte Carlo simulations of 55 random field positions (equal to the 
number of unidentified SMGs), and employ the mean and standard 
deviation of the 50 simulations in each redshift bin for our statis- 
tical analyses. As discussed in Sj2. 1 l our photometric source extrac- 
tion procedure included manually examining the regions around the 
SMGs and adding to the catalogue potential sources which may 
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have been missed by the automated procedure. To remove any bias 
and ensure a fair comparison between the SMGs and random posi- 
tions, we exclude these additional sources from this analysis. 

In the lower panel of Fig. [5] we show the difference between 
the redshift distributions of the field and SMGs without robustly 
identified counterparts. Compared to the field there is an excess 
of 26 ± 12 (2.2<t) at z > 1 around unidentified SMGs. There are 
positive excesses of galaxies in the z = 2-3 (14±8; 1.8cr) and z = 
4-5 (4 ± 2; 1.6a) bins around unidentified SMGs. This suggests 
that the peak of the redshift distribution of the population of radio, 
24-^.m and IRAC unidentified SMGs is at z = 2-3. 

To crudely compare the redshift distributions of the identified 
and unidentified SMGs we also plot the difference between red- 
shifts of sources in the submillimetre error circles of the identified 
SMGs to the field (scaled to the value in the z = 2-3 bin of the 
unidentified SMGs). We conclude that the redshift distribution of 
unidentified SMGs is broadly similar to that of robustly identified 
SMGs. 

To provide a more reliable estimate of the average redshift 
of the unidentified SMGs we evenly distribute the excess galax- 
ies in the SMG error circles in each redshift bin. We verify that 
this method is valid by using it to calculate the average redshift 
for identified SMGs, which yields z = 2.2 ± 0.1, in agreement 
with that derived using the robust counterparts alone ( i]4.1t . For 
unidentified SMGs we derive an average redshift of z = 2.5 ± 0.3. 
This suggests that unidentified SMGs may lie at marginally higher 
redshifts than the identified sample, although we stress that the 
difference is not statistically significant. This conclusion is con- 
sistent with the predicted redshift distribution of radio undetected 
SMGs from the semi-analytic ACDM GALFORM model, which 
predicts they should li e at z ~ 2.2, similar to the observed SMGs 
dSwinbank et al J2008t) . 

There are 14 ± 8 more galaxies at z — 2-3 in the error circles 
of unidentified SMGs than expected from comparing to the field 
population. We showed in £|3.3I that there are nine tentative SMG 
counterparts with z = 2-3 and two with z — 4-5. Thus poten- 
tially half of the excess seen around the unidentified SMGs could 
be attributed to tentative counterparts. Indeed, if tentative SMGs 
are removed from this analysis a 0.6ct (3 ± 5) excess of z — 2-3 
galaxies and a 1.3<r (2.4 ± 1.9) excess of z — 4-5 galaxies around 
SMGs remains. 

Finally, we use statistical arguments to estimate how many 
SMGs are still unaccounted for. There are 55 out of the sample 
of 126 LESS SMGs wit hout robust radio, 24 y.m or IRAC counter- 
parts jBiggs et"aT1l2010l) . Due to the signal-to-noise ratio limit on 
the submillimetre catalogue (S/N > 3.7a) five of the 126 SMGs 
are expected to be false detections JWeifl et alj2009h ; an additional 
1-2 a re expected to hav e counterparts outside of the search radii 
used jBiggs et al1l2010h . This leaves 48-49 SMGs with currently 
unidentified counterparts that are expected to lie within our search 
area. We then calculate the total excess of galaxies in unidenti- 
fied error circles over the field. Due to clustering, an error circle 
can contain more than one galaxy associated with the SMG. We 
then compare the sources around identified SMGs with the number 
of identified counterparts to determine this "overcounting factor": 
~ 1.2 x. We scale the difference between the field and unidenti- 
fied SMG regions by this factor to estimate that there are 21 ± 19 
LESS SMGs (17 ± 15% of the total) which are still unaccounted 
for. These have no robust radio, 24 \im and IRAC counterparts and 
have mid-infrared fluxes below the limits of our imaging. 

These counterparts could lie at z ~ 1-3 and be fainter than 
Mh < —23 (Fig. 1 10b . However, the specific star-formation rates 



of such sources would be > 10 yr _ , corresponding to life 
times of < lOMyr. The corresponding duty cycle for such short- 
lived sources means that to detect ~ 20 sources in our 0.5 deg 2 
survey, we require a parent population with a space density of 
> 0.02 Mpc , which we consider unlikely. Alternatively, if they 
have restframe near-infrared luminosities similar to the identified 
SMG population, then Fig.[l0]suggests that they must lie at z > 3. 
If correct we should add these sources to the SMGs identified 
at z > 3. We have identified 10 SMGs at z > 3, as well as 
4 ± 2 which have been statistically identified in our IRAC sam- 
ple. To these we add the 21 ± 19 SMGs which are unaccounted 
for in our statistical analysis, to derive a total of 35 ± 19 SMGs 
(28 ± 15% of the whole population) at z ^ 3 in our survey. We 
conclude that ~ 30%, and at most ~ 45% of the SMG popula- 
tion could reside at z > 3. This corresponds to a volume density 
of 2.8 x 10 -6 Mpc -3 (assuming they span the range z — 3-7, 
or 80% higher if they only span z = 3-5). For comparison, the 
equivalent volume density of z = 2-3 SMGs, including identified 
counterparts and the 14 that are statistically identified in this red- 
shift range, is 1.2 x 10~ 5 Mpc -3 , signifying strong evolution in the 
abundance of SMGs from z > 3 to z ~ 2.5. 

We have statistically identified the redshifts of ~ 60% of 
the unidentified SMGs, and shown that the remainder likely lie at 
z > 3. In the right-hand panel of Fig. [3] we combine the redshift 
distributions of the identified and unidentified SMGs to provide the 
probable redshift distribution of the entire Saroum > 4 mJy SMG 
population. We conclude that the most likely median redshift for 
the Ss7o^m > 4 mJy SMG population is z = 2.5 ± 0.6. 



4.3 Simple redshift estimators for SMGs 

Previous studies h ave investigated and used optical (BzK ), ul- 
traviolet (BX/B M; Isteidel et al l 12001: IChapman et all 120051) and 
IRAC colours jYun et alj|2008l : iHainline et alll2009l) and radio- 
to- submillimetre flux ra tios (e.g. ICarilli & Yunll 19991 : llvison et al.1 
120071 : iBiggs et alj|2010h as simple estimators of the redshifts of 
SMGs. Here we use our 17-band photometric redshifts to in- 
vestigate the reliability of such the BzK colours and radio-to- 
submillimetre fluxes as redshift estimators and derive a simple 

IRAC colour indicator of redshift. 

In Fig. [6] we show the BzK colour-colour plot jDaddi et al.1 

2004), which is designed to identify galaxies at 1.4 < z < 2.5, 
for LESS SMG counterparts. We have distinguished between coun- 
terparts with photometric redshifts above and below z = 1.4 and 
find that all the SMGs with 2 p i 10 t < 1-4 lie in the expected region 
of colour-colour space. However, whilst SMG counterparts with 
Zphot 1.4 typically have the colours of high-redshift star-forming 
galaxies, this population does scatter into the low-redshift region. 
Two galaxies with z p hot > 1.4, B — z ~ 1 and z — K ~ —0.1 
lie near the separation between 2 < 1.4 galaxies and stars and 
are both X-ray luminous, bright unresolved sources which may be 
submillimetre-bright quasars (see Appendix [A}. We conclude that 
the BzK analysis of SMG counterparts can select clean but in- 
complete samples of z > 1.4 SMGs, and that samples selected to 
have z < 1.4 will contain ~ 45% contamination from galaxies 
at higher redshift. We compare the observed BzK colours with a 
redshift track of the average SMG SED ( jgjD and note that the 
median SED of SMGs has a redder restframe (U — 2) colour (cor- 
responding to observed (2 — K) at 2 ~ 1.4) than used to define 
the selection areas for 2 > 1.4 galaxies and that the SMGs at the 
highest redshifts may fall in the passive BzK or 2 < 1.4 region. 
In Fig. [7] we plot the photometric redshift against 



© 2010 RAS, MNRAS 000, [711231 



12 J. L. Wardlow et al. 




Averoge SMG track < 



1 



2 3 4 

B-z 



Figure 6. (B — z) versus (z — K) colour-colour plot of LESS SMG 
counterparts. The selection regions for star-forming and passive z > 1.4 
BzK galaxies (sBzK a nd pBzK respectively), z < 1.4 galaxies, and 
stars toaddietal]|2004l) are shown, and we distinguish between SMG 
counterparts with z p hot ^ L4 and z p i lot < 1.4. The photometric red- 
shifts typically agree with the BzK colours and most of the SMG coun- 
terparts have BzK colours of z > 1.4 star-forming galaxies and none 
have colours of star s or z > 1.4 passive galaxies (similar to the result of 
iBertoldi et al]|2007t) . All of the counterparts with sBzK colours are found 
to have z p hot ^ 1.4, but galaxies with BzK colours suggesting z < 1.4 
have a ~ 45% contamination from counterparts with z p ^ ot ^1.4. We also 
show the redshift track of the average SMG SED W.4t from z = 0-4 and 
the reddening vector for Ay = 1 mag. The scatter in the photometry of 
individual SMGs compared to the redshift track of the average SMG SED 
suggests that the SMGs have a range in optical SEDs. Open symbols show 
galaxies which lie in halos of bright stars in the z-band, in these cases, for 
the purpose of this plot only, the z-band magnitude is extrapolated from the 
SED fit and the measured /-band magnitude. 



S870um/Si . 4GHz for the LESS SMGs jWeiB et alj 120091 : 
iBiggs et al.l l2010h. the tracks of Arp 220 and M82 (based on 



the SEDs of lSilvaetal1ll998l). We also show the ACDM gal- 
form predictions teaugh et ai] [2005) and the ICarilli & Yunl 
(2000) relationship. The wide range in S's70m™/'S'i.4GHz at a fixed 
redshift limits the usefulness of SWonm/Si^GHz as a re dshift 
indicator for SMGs and indicates that SMGs have a variety of 
submillimetre -to-radio flux ratios , suggesting a range in dust 
temperatures ^Chapman et al.l2005l : lciements et alj2008l) . We also 
note that the majority of SMGs lie above the redshift track of M82, 
suggesting higher submillimetre-to-radio flux ratios (potentially 
due to the presence of more cold dust). LESS 20 has z p hot ~ 2.8 
and S , 870/jm/S'i.4GHz ~ 1.7, which is significantly lower than 
expected from its redshift, indicating that it is most likely a 
radio-bright AGN, so we remove it from our subsequent analyses 
of far-infrared luminosities, star-formation rates, and characteristic 
dust temperature (i ]4.6l >. 

Studies of mid-infrared spectra of SMGs have shown that they 
are similar to M8 2 with an additional power-law contribution from 
AGN emission |M enendez-Delm estre et al.l2007l.l2009MPope et al.l 
l2008l : ICoppineta l. 2010b). In Fig. [7] we plot photometric redshift 
against iSWo/tm/SW^ra, wmcn shows that the mid-infrared to sub- 



millimetre flux ratios of SMGs are similar to Arp 220 and those 
derived in the ACDM GALFORM model, but are poorly represented 
by M82. This suggests that although SMGs have mid-infrared spec- 
tra similar to M82, the mid-infrared continuum emission is fainter 
compared to the far-infrared emission and is more comparable to 
that of Arp 220. We note that although SWo^m /SW^m for SMGs 
varies with redshift in a manner comparable to Arp 220, the scat- 
ter and the effect of PAH and silicate features passing through the 
24 \i m filter makes thi s measurement unsuitable for redshift deriva- 
tion JPope et al.ll2006l) . 

W e expand on the work o flYun et al.1 J2008h and lHainline etail 
d2009h and propose a new redshift estimator for SMGs, which is 
based on the IRAC 8 and 3.6 \im fluxes and exhibits less scatter 
than the commonly-employed radio-to-submillimetre flux ratio. In 
Fig. [8] we plot this ratio against redshift for the LESS SMGs and 
using the ROBUSTXIN EFIT procedure from the IDL Astronomy Li- 
brary ([Landsman 1993) we fit an outlier-resistant linear relationship 
to SMGs with Zphot < 4, which yields: 



= (2.1 ± 0.1) + (1.9 ± 0.2)log 10 (5 8 /S 3 . 6 ) 



(1) 



We exclude SMGs with z p i lot ^ 4 from the fit because at high red- 
shifts the 1.6 \im stellar peak passes through the 8 urn filter making 
this redshift estimator unreliable. For SMGs with z < 4 the la 
dispersion in redshift estimated using Eqn.Q]is a z — 0.44 and we 
find that ~ 90% of SMGs at z > 2 have S 8 /S 3 . 6 > 1. 

4.4 Typical SMG SEDs 

To investigate the SED of a typical SMG we show in Fig. [9] the 
SMG photometry in the rest-frame and normalised in the 17-band. 
We also calculate the expected fluxes expected in each of the 17 
photometric filters used throughout this paper as observed at z = 
2.2 and determine the median flux in each. 

We begin by noting that the data show evidence for a break at 
~ 3500-4000A suggestive of a Balmer or 4000A break. Closer 
inspection hints at it being a Balmer break indicating that the 
blue rest-frame light is dominated by stars older than 20 Myr and 
younger than ~ 2Gyr. Then, as in ij3.ll we use HYPERZ to fit 
this photometry, with redshift fixed at z — 2.2, and show the best- 
fit templates for both the Burst and Im star-formation histories on 
Fig. [9] Comparing the \ 2 for these two models in the same man- 
ner as i]3.5l we find that we cannot accurately distinguish between 
different star-formation histories (and hence ages or light-to-mass 
ratios). The best-fit Burst model has an age of just 33Myrs, an 
Ay = 1.7 and a resulting light-to-mass ratio of Lh/M* ~ 24, 
in contrast the Im template yields an age of 3.5Gyrs, Ay = 1.1 
and an Lh /A/* ~ 6. The reddening derived from these two tem- 
plate fits are agreement with the median of the individual SED fits 
(Ay = 1.5 ±0.1; CHi while the Lh/M* span a range of 4x. 

We also estimate the extinction in LESS SMGs by compar- 
ing the star-formation rate (SFR) derived from the rest- frame far- 
ultrav iolet luminosity (median SFRuv = 2 M0vr~ 1 ; lKennicuttl 
1 19981) with the SFR derived from the far-infrared luminosity (me- 
dian SFRfir = 1400 M yr _1 ; H6I >. Comparing the two values 
yields Ay = 2.6 ± 0.2, corresponding to reddening ~ 4 times 
higher than the SED fit and indicating that the majority of the star 
formatio n within SMGs oc curs in totally obscured regions. As dis- 
cussed in lKennicuttl (1998) the conversion from far-ultraviolet lu- 
minosity to SFRuv assumes that the star-formation rate has been 
constant for > 10 8 years. SMGs are likely to be shorter bursts of 
activity and therefore for a fixed SFR they will brighter at ultravio- 
let wavelengths and likely have higher Ay than estimated above. 
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Figure 7. The variati on of submillimet re to radio (left) and submillimetre to mid-i nfrared (right) flux ratios with redshift, compared with Arp 220 and M82 
(based on the SEDs o flSilva et aTJ l998) and SMGs in the ACDM GAL FORM model jBaugh et al.l20bH) with SssOfim > 3 mjy; we also show the relationship 
between redshift and radio-to-submillimetre spectral index derived bv lCarilli & Ymil | 2000|) in the left-hand panel. The SMGs show two orders of magnitude 
dispersion in both 5870nm/Si.4GHz and S870/jm /S24um ■ The model track of 5"870a™ /'-'1.4GHz for M82 lies below the majority of the SMGs, while that 
of Arp 220 more closely follows the SMGs suggesting that they typically have a characteristic dust temperature which is cooler than that in M82 and more 
like t hat in Arp 220. Similarly, although studies have found that mid-infrared spectral properties of SMGs are similar to M82 ( Menendez-Delmestre et al. 
|2009|) we find that M82 does not describe the submillimetre to mid-infrared continuum flux ratios well, and that Arp 220 fits better to this data. LESS 20 
(labelled) is significantly brighter at 1 .4 GHz than expected from its submillimetre flux and redshift and is most likely a radio-bright AGN. 




Figure 8. The correlation between redshift and the ratio of 8 |xm to 3.6 \im 
flux for the LESS SMGs. This shows a trend and so we plot a linear fit to 
the SMGs with 2 < 4, which yields 2 = 2.1 + 1.91og 10 (S l 8/>S3.6). with 
a 1<t dispersion in redshift of tr z = 0.44. This relation may be useful as a 
crude redshift indicator for SMGs as we note that ~ 90% of all SMGs with 
2 > 2 have 58/53.6 > 1, while similarly ~ 90% of all SMGs with 2 < 2 
have S$/S 3 .G < 1. 

4.5 Stellar masses 

We use HYPERZ to estimate the rest-frame //-band absolute mag- 
nitudes (Mh) from our SED fits and find the median Mb = 
—24.1 ± 0.1, with an interquartile range of —24.7 to —23.6. In 
Fig. [TO] we plot Mh against photometric redshift for the LESS 
SMG counterparts. There is the suggestion of a weak trend of Mh 
with redshift. However, as the plotted detection limit shows this 
is most likely a selection effect, with the higher redshift galaxies 
needing to be more luminous to be detected. The average z = 2-3 
SMG has Mh — —24.4 and would be detected in our survey out 
to z ~ 3^1, and therefore, any incompleteness in our SMG sample 
due to the IRAC flux limits is only significant at z > 3. 

Mh is used to estimate the stellar mass of galaxies because it 



Figure 9. The photometry of SMG counterparts shifted to the rest frame 
and normalised to the H band (1.6um), Redshifting this data to 2 = 2.2 
we calculate the apparent fluxes in the 17 photometric filters considered 
throughout this paper and use HYPERZ to fit galaxy templates at this red- 
shift. The median photometric points are shown and the resulting best HY- 
PERZ Burst and Im template fits are displayed. The best fit HYPERZ tem- 
plates have: Burst, Ay = 1.7, an age of 33Myrs and a Lh/M* ~ 24; 
Im, Ay = 1.1, an age of 3.5 Gyrs and an Lh /M* ~ 6. It is clear that it is 
not possible to distinguish between these two very different star-formation 
histories and hence there is a factor of ~ 5 uncertainty in the resulting 
masses. The MUSYC U38 filter has an ^ 50% contribution from limiting 
magnitudes and is excluded from the fit. 

is less influenced by young stars than optical bands and is relatively 
unaffected by dust. As discussed in £13.51 the uncertainties in the 
derived spectral types and ages result in an estimated factor of ~ 5 
uncertainty in assumed mass-to-light ratios and thus stellar masses 
derived from Mh - Therefore, we only consider the stellar masses 
of t he LESS SMGs statisti cally. 

lHainline et all d2010h estimated //-band ma ss-to-light ratios 
for SMGs with Burst and Im templates, based on a lChabrieil d2003h 
IMF . We use the average of their values converted to a Salpeter 
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Figure 10. A plot of photometric redshift against rest-frame //-band abso- 
lute magnitude for LESS SMGs, and the approximate correspondence with 
stellar mass (as described in the text). The median Mjj is —24.1 ±0.1 
with an interquartile range of —24.7 to —23.6 which corresponds to a 
median stellar mass of ~ 9.2 X 10 10 Mq and interquartile range of 
(4.7-14) X 10 10 M (Ti in good agreement with the median Mh from 
lHainline et allfeOlOl) . We also highlight SMGs with evidence for AGN ac- 
tivity from X-ray detections or 8 urn excesses, which appear brighter than 
the average SMG. We also show the 1-cr distribution of absolute //-band 

3 mjy that are brighter than our 
2) from 



magnitudes of SMGs with S850^« 

flux limit a t 4.5 txm (approximately the rest-frame //-band at z 
GALFORM iBaugh et alj2005tlLacev et alj2008r . lSwinbank et al.l2008l) and 
note that Mh is under-predicted in the model. Errors in Mh are dominated 
by the error in the photometric redshift; we calculate the error in Mh for 
SMGs with median redshift error by re-running the HYPERZ with the red- 
shift forced to the extremes of the error range; the corresponding errors are 
shown in the top left-hand corner of the plot and can be scaled with the 
error in redshift. The dotted line illustrates the trend in Mh with redshift 
resulting from the flux limited nature of our survey. 



IMF (with a lower mass limit of 0. 1 Mq and an upper mass 
cutoff of 100 Mq) for our stellar mass calculation: Lh/M* = 
3.8 LqMq 1 . We estimate that the median stellar mass of the 
SMGs in our sample is M* = (9.2 ± 0.9) x 10 10 Mq and the in- 
terquartile range is (4.7-14) x 10 10 Mq. The quoted errors do not 
include the systematic uncertainty from the star-formation histories 
and mass-to-light ratios, which adds a factor of ~ 5 uncertainty to 
the values ( ^3.51 Fig. |9). We also caution that the choice of IMF 
coupled with the assumption that all the light is from the current 
burst can affect the derived stellar masses by an additional factor 
of ~ 2. Finally, we note that on average we observe the SMGs ap- 
proximately halfw ay through the burst and typical SMG gas masses 
i lGreve et al]|2005h suggest an additional ~3x 10 10 M Q could be 
added by the end of the burst. 

We find that galaxies with evidence for AGN activity from an 
8 [im excess or X-ray emission have median Mh = —24.6 ± 0.3, 
compared to Mh = —24.1 ± 0.1 for the remainder of the SMGs. 
The two SMGs with the brightest Mh are the two submillimetre 
bright quasars (LESS 66 and LESS 96; Appendix [At in which the 
observed emis sion is expected to b e dominated by the AGN rather 
than starlight dHaimineetall koiO). If these are excluded the me- 
dian M H of SMGs containing AGN is M H = -24.5 ± 0.3. 

The median stellar mass for SMGs in SHADES Lockman 
Hole was claimed to be M* = (6.3+i g) x lO n M bv lDve et al.1 



(2008). This is a factor o f ~ 7 higher than our estimate for LESS 
SMGs. iDve et all d2008f) use nine-band photometry for their pho- 
tometric redshift determination and claim to also be able to dis- 
entangle the star-formation histories of the SMGs with sufficient 
accuracy to identify a significant mass of old stars which underlies 
the current burst. This leads to a higher effective mass-to-light ra- 
tio and correspondingly higher stellar masses. In contrast, as dis- 
cussed earlier ( i]3.5t . we do not believe that with existing data 
it is possible to untangle the influences of the potentially com- 
plex star-formation histories and dust distributions on the SEDs of 
SMGs. Hence, we do not believe that there is any observational 
evidence for signi ficant old s tellar populati ons in these galaxies , 
as required by the lDve etal] d2008l) results. lHainline et al.l feoifj) 
have used optical and IRA C photometric data to calculated an av- 
erage stellar mass for the IChapman et all (2005) SMGs and they 
find M* = (1.4 ± 0.3) x lO n M (converted to Salp eter IMF) 
compa rable to our survey and a factor of ~ 5 lower than lDve et al.l 
d2008l) . 

In Fig. [10] we also show the absol ute //-band magni tudes of 
SMGs in the ACDM GALFORM model dBaugh et al]|2005h . which 
assumes a top-heavy IMF with slope x — 0. We consider only 
galaxies with Ssso^m 3 mjy and fluxes in the IRAC 4.5 [im fil- 
ter brighter than our dete ction limit (4.5 yim corr esponds to the rest- 
frame //-band at z ~ 2). ISwinbank et al! I d2008l) showed that GAL- 
FORM predicts rest-frame absolute A"-band luminosities of SMGs 
which are a factor of ten lower than observed. This arises primar- 
ily due to an order of magnitude lo wer stellar masses than implied 
by observations for SMGs (see also lLacev et alj |2010). As Fig. 1 101 
shows the predicted rest-frame //-band magnitudes of the model 
SMGs are also a factor of ten lower than our observations. In- 
deed, if SMGs for med stars following the prescriptions used in the 
iBaugh et al.l d2005l) . then few of the SMGs above a redshift of z ~2 
would have be detected. 



4.6 Dust temperatures, far-infrared luminosities and 
star-formation 

In order to further investigate the intrinsic properties of the LESS 
SMGs we next use our photometric redshifts and the observed ra- 
dio and submillimetre fluxes to derive the characteristic dust tem- 
peratures (Td), far-infrared luminosities (8-1000 [im; Lfir) and 
star- formation rates. 

Blai n et al. I d2002l) showed that the submillimetre-to-radio flux 
ratio in SMGs is mainl y influenced by redsh ift and the characteris- 
tic dust temperature. IChapman et al.l q2005) assumed a dust emis- 
sivity, P = 1.5, and the z = far-infrared-radio correlation, to 
determine empirically that for their sample of SMGs: 



T D = 



6.25(1 + 2) 



(S'850Aim/5'l.4GHz) 



(2) 



We note that the most reliable method of calculating Td is to fit 
template SEDs to multiple far-infrared and submillimetre photo- 
metric points, but for simplicity and due to the absence of pub- 
lished deep far-infrared photometry we use Eqn.[2]to calculate Td 
of LESS SMGs (although we next use shallow far-infrared obser- 
vations to confirm the validity of this assumption). 

We also use the infrared-radio correlation dHelouetalJI 19851: 
ICondon|[l992l) , 



qFIR = log x 



Li 



,3.75 x 10 12 W 



logx 



f- 

Vw 



WHz- 



(3) 
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Figure 11. The characteristic dust temperature (To) versus far-infrared lu- 
minosity (Lpjr) for our SMGs. The SMGs are colour coded on the basis 
of their photometric redshifts and as expected the most luminous galaxies 
are the hottest, and are also tend to be those at the highest redshifts. This 
trend is driven in part by the radio luminosities of the SMGs (which lack the 
positive K correction of the submillimetre waveband) hence why there is a 
correlation between Lpm and z, but not between S , 870/jm and z (Fig. |4}. 
The regions above the dashed line and below the dotted line are illustrative 
of the regions excluded by our submillimetre and radio detection limits re- 
spectively. The dashed line, which roughly demarcates the upper envelope 
of the data, represents the derived temperature of galaxies at z = 2 with 
various radio fluxes and SsTa^m = 4.2 mJy. The dotted line, which simi- 
larly demarcates the lower envelope, is derived for submillimetre-luminous 
(Sg70/j,m = 16 mJy) sources with radio flux equal to our detection limit 
(3<t = 19.5pJy) at redshifts of z = 1.4 — 4. This reflects both the strong 
cut-off in the submillimetre luminosity function at high luminosities and 
the fact that our radio data is only just deep enough to detect counterparts to 
the majority of SMGs. We conclude that the apparent correlation between 
Tp and L fir is in part caused by selection bias. We note that the OFRGs 
(Chapman et al. 20043: ICasev et al.lEo 09: Magn elli et ai]|2010l) . which are 
detected in the radio but not the submillimetre lie above the upper dashed 
line. 



Figure 12. Far-infrared luminosity functions of the radio-detected LESS 
SMGs with z = 1-2 and z = 2-3. Evolution is evident in the luminosity 
function of the LESS SMGs between the two re dshift bins, and from th e 
24- urn selected low redshift comparison sample I Rodighi ero et al ] l201Ch . 
The SMGs have higher luminosities than the z < 0.3 24-um galaxies, 
and the z = 2-3 SMGs having higher luminosi ties and <E>* . We also show 
the z = 2-3 luminosity function of IChapman et alj (2005) SMGs (offset 
slightly in log 10 Lpjn for clarity) for comparison. The L ESS z = 2-3 
sampl e has a systematically lower lumino sity density t han|C hapman et al. 
12005) SMGs in the same redshift range. IWeiB et alj (2003) " showed that 
the ECDFS is underdense at submillimetre wavelengths and by scaling the 
LESS luminosity fu nction such tha t the SMG number density matches that 



of SHADES survey I Coppin et al 



between LESS and Chap man et al 



2006) we show that the disparity in $ 



i2005h is likely due to the relative den- 



sity of SMGs in the two surveys (the scaled LESS luminosity function is 
offset slightly in \og 10 LpiR SMGs for clarity). We calculate the maxi- 
mum contribution from unidentified SMGs, by assigning them the redshift 
distribution that we measure in ij4.2| and radio fluxes equal to our detection 
limit. Including the contribution from unidentified SMGs the total maxi- 
mum $ in each luminosity bin is represented by an arrow (offset slightly in 
log 10 Lj?/ij for z = 2-3 SMGs for clarity). 



with radio spectral index a — 0.8 (where Su oc u~ a ) and 
<1fir = 2.64 (Bell 2003, for star-forming galaxies), to calculate 
far-infrared luminosities of the LE SS SMGs from their radio fluxes, 
as done by Chapman et al. ( 2005). Although this approach was re- 
cently verified by Magnelli et alj (20 id) who used Herschel data 
to show that the local far-infrared radio correlation is consistent 
with SMGs, we caution that there may be a factor of ~ 2 uncer- 
tainty in the derived luminosit ies due to possible evolution in the 
far-infrared-radio correlation (ivisonetal]|2010al) and hence the 
appropriate value of qfir- 

In Fig. QT] we plot the far-infrared luminosity against To 
for the LESS SMGs and optically faint radio galaxies (OFR Gs; 
IChapman et al.l2004al : lc"asev et all 20091 ; iMagnelli et alfcOld . The 
LESS SMGs have a median T D = 35.9 ± 1.4 K, with an interquar- 
tile range of 28.5-43.3 K, and L FIR = (8.2 ± 1.2) x 10 12 Lq, 
with an interquartile r ange of (3.0-13) x 10 12 L^, comparable to 
previous surveys (e.g. Chapman et al. 2005; Magnelli et alj|201ol) . 

To check this result we also employ the 250, 350, and 500 \im 
Balloon-borne Large Aperture Submill imeter Telescope (BLAST) 
maps of the ECDFS (Devlin et alj2009h . We can stack the emission 
in these maps at the positions of the LESS SMG counterparts and 
fit the stacked fluxes with a modified black body with ft — 1.5 at 
z — 2.2 and correct the luminosity of the fit ted black body to to - 
tal infrared luminosity, 8-1000 fim, based on llvison et alj (2010ch . 
From this calculation the typical characteristic dust temperature of 



the LESS SMGs is Td = 33. 6± 1.1, and the typical far-infrared lu- 
minosity is Lfir = (7.&t\' B ) x 10 12 Lq. These values are in good 
agreement with those derived above from the local far-infrared- 
radio correlation. 

We find that the highest redshift galaxies also have the highest 
luminosities due to a combination of the radio K-correction (pre- 
venting the detection of low-luminosity galaxies at high redshifts) 
and luminosity evolution (see Fig.ll2t. There is an apparent trend 
between the To and Lfir but this is likely at least partially a se- 
lection effect, although we note that locall y IRAS galaxies exhibit 
a tight correlation between Td and Lir (Chapman et al. 2003b; 
Chapi n et alj2009l) . To illustrate the selection effects we also show 
in Fig.ll HOFRGs ( Chapman et al. 2004a), which are detected at ra- 
dio but not submillimetre wavelengths and have radio luminosities 
similar to SMGs, but contain warmer dust (Td ~ 45 K; Ca sey et all 
120091 ; iMagnelli et aljfeoiol) . Our 870-nm detection limit misses 
warmer and lower luminosity galaxies from the sample and the 
radio detection limit excludes the colder luminous galaxies (e.g. 
IChapman et"ai]|2005l) . 

In Fig. [12] and Table [3] we present the far-infrared luminosity 
functions of the radio-detected LESS SMGs wit h z = 1-2 and 
2 = 2-3, compared to the z = 2-3 result from Chapm an et al.l 
(2005). We calculate the LESS SMG luminosity function with an 
accessible volume technique where: 
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Table 3. Far-infrared luminosity function for radio-detected LESS SMGs 



z = 

logioI/FIfl 
(£©) 


1-2 SMGs 

logio $ 
(Mpc- 3 dcx _1 ) 


3 = 

logioI/F/fl 
(%) 


2-3 SMGs 

logio $ 
(Mpc~ 3 dex _1 ) 


12.0 
12.5 
13.0 
13.5 


r ct+0.2 
0.3 

-5.3 ±0.1 

-5 3+ - 1 
O.J_ 2 

< -6.2 


12.5 
12.9 
13.3 
13.7 


°' 1 -0.6 

-5 1+° 1 
°' 1 -0.2 

c q+0.2 

-6 2 +a7 



(4) 



LESS SMGs f 
Chapman et al 2005 I 

LESS scaled to SHADES SMG density ♦ < 

Hopkins & Beacom 2006 

Evolution from z = 




which accounts for the flux limited nature of our survey. $(L) AL 
is the number density of sources with luminosities between L and 
L + AL, and Vi is the comoving volume within which the ith 
source can be detected in the luminosity bin under consideration. 
Since we derive the far-infrared luminosity from the radio flux, Vi 
is calculated using the radio luminosity. Error bars are calculated 
by bootstrapping and account for the redshift, luminosity and bin- 
ning errors. We use the same meth od to calculate the luminosity 
function for Chapman et alj d2005l) SMGs based on the redshifts 
and radio fluxes listed in that paper. By assuming that unidenti- 
fied SMGs have radio fluxes equal to our detection limit and the 
redshift distribution that we measure in t\4.2\ we also calculate the 
maximum contribution of unidentified SMGs to the far-infrared lu- 
minosity functions. 

We observe strong evolution in the far-infrared luminosity 
function: the z — 2-3 SMGs are more luminous and have 
higher space densities than the z — 1-2 SMGs, which in turn 
are more luminous th a n the z < 0.3 24 itm-selected galaxies 
from [R odighiero' et alj d2010h (see also the Herschel sample of 
IVaccari et alj|2010t) . LESS SMGs at z = 2-3 have $* ~ 60% 
and L* ~ 2x larger than those at z — 1-2. 

The z = 2-3 LESS S MGs have systematically lower <J>* than 
the lChapman et alj J2005I) SMGs i n the same redshif t range. This 
may be due to cosmic variance since IWeiBetalJfc009h showed that 
the ECDFS is a factor of ~ 2 underdense compared to other large 
submillimetre surveys at flux densities > 3 mJy. By rescaling the 
LESS luminosity function so that the 870 p.m number coun ts agree 
with those of the SHAD ES survey dCoppin et alj feoOfj) . which 
should be similar to that of Chapman et al. (2005) since both cov- 
ered multiple fields, Fig. [T2] shows that the low surface density of 
SMGs in the ECDFS is most likely the cause of the disparity in $* . 

In Fig. [T3] we show the evolution of the star-formation rate 
density (SFRD) of the radio-detected LESS SMGs. We use the 
same accessible volume technique as in our luminosity function 
calculations to account for the flux limited nature of the survey. 
Error bars are calculated from bootstrapping and include the uncer- 
tainties in binning, redshifts and SFRs. Since the SFRs are based 
upon radio fluxes we exclude the suspected radio-bright AGN 
LESS 20 from this analysis. 

We do not know the individual redshifts, infrared luminosities 
or SFRs for 45% of the LESS SMGs because they do not have ro- 
bustly identified optical counterparts. In Fig.[T3]we account for this 
population by assigning them the redshift distribution that we mea- 
sure in S |4.2b and assuming radio fluxes equal to our detection limit. 
The calculated SFRD of the unidentified SMGs from this analysis 
is an upper limit since the actual radio fluxes will typically be lower 
than the detection limit. In Fig.[T3]we indicate the maximum con- 
tribution to the SFRD of unidentified SMGs in each redshift bin. 



Figure 13. Evolution of the SF RD for the radio-detected LESS SMGs com- 
pared to lChapman et alj (2005). Arrows to the right of each LESS redshift 
bin indicate the maximum additional contribution from unidentified SMGs 
and th e open symbol represents unidentified SMGs from Chapman et al. 
t2005t) . We a lso show the modified Salp eter A IMF fit to the SFRD com- 
pilation from IHopkins & Beacom (2006) a nd a line show ing the evolution 
from IRAS ULIRGs at z = lEfbaz & Cesarsky||2003h to LESS SMGs 
at z = 2.3. The LESS SMG activity peaks at z ~ 2 - similar to that 
found by pr evious studies of s tar-forming galaxies and the peak activity 
of quasars (Hopkins et al. 2007j). The contribution from SMGs to the to- 
tal SFRD also peaks at z ~ 2 where they are responsible for ~ 10% of 



the IHopkins & Be acom 



20061) SFRD. Th e ECDFS is underdense at sub 



millimetre wavelengths WeiB et al. 2009) so similarly to Fig. [12] we also 
scale the SFRD of the LESS SMGs such that the number counts match the 
SHADES survey allowing a closer comparison to Ch apman et alj ^2005). 



The SFRD of the LESS SMGs appears to peak at z ~ 2, simi- 
lar to lChapman et al.1 J2005I). The LESS SMG s have a lower SFRD 
than the SMGs from Chapman et alj (2005) but we note that the 
lower number density of SMGs in the ECDFS is sufficient to ac- 
count for this effect. This corresponds to t he peak of quasar activity 
at z — 2.15 ± 0.05 (Hopkins et al. 2007). The fractional contribu- 
tion of LESS SMGs to the SFRD of the Universe also peaks at 
z ~ 2 w here they are responsible fo r ~ 10% of the SFRD as esti- 
mated bv lHopkins & Beacoml (2006J) from a compilation of surveys 
that does not include any submillimetre surveys. We stress that this 
only includes SMGs with SsTOfim 4 mJy. Assuming that fainter 
sources have the same redshift distribution then the contribution of 
SMG s with SWrjAmi > lmJy is ~ 100% of the (Hopki ns et~ai1 
120071) value. Thus, SMGs contribute ~ 50% of the total SFRD of 
the Univers e at z ~ 2. 

We use iKennicuttl fl998l) , which assumes a Salpeter IMF with 
upper and lower mass limits of 0. 1 and 100 Mq respectively, to cal- 
culate the SFRs of the LESS SMGs from their inferred far-infrared 
luminosities. The median SFR is 1100 ± 200 MQyr - and the in- 
terquartile range is 300-1900 MQyr -1 . The median specific star- 
formation rate sSFR = SFR/M* = (1.2 ±0.1) x lO^V" 1 and 
the interquartile range is (0.6-1.8)xl0 -8 yr~ 1 - Although again, we 
caution that due to the uncertainties in the stellar mass estimates 
there is an additional factor of ~ 5 uncertainty in these values (see 
^3 .51 for a full discussion). The median formation timescale of the 
LESS SMGs is thus ~ 100 Myr and it is feasible that all of the stel- 
lar mass we see could be formed in a single burst. In Fig.[l4]we plot 
the trend of sSFR against the redshift of the LESS SMGs. We note 
that the apparent lack of galaxies with low sSFR at high redshifts is 
a selection effect due to the requirement for a radio counterpart and 
that galaxies with sSFR > 10 -7 yr -1 are rare due the brevity of the 
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Figure 14. A plot of specific star-formation rate (sSFR) versus redshift 
for the LESS SMGs. Galaxies are colour-coded by mass and we show the 
median error bar in the top-left hand corner and note that similarly to Fig. 1 101 
the error in sSFR is correlated with that in redshift. The arrow represents the 
gradient of the trend in sSFR with redshift for IRAC-selected galaxies with 
log 10 (M*) = 10.3-10.8 Mq, offset in sSFR by two orders of magnitude 
for the purpose of display. We note that due to the requirement for radio 
counterparts no SMGs are detected in the high redshift and low sSFR region 
of this plot. Similarly, the short lifetime of SMGs with sSFR > 10 _7 yr _1 
and the limited volume of our survey means that few SMGs with very high 
sSFRs are detected. However, the dearth of SMGs at 2 < 1.5 with sSFR 
> 10" 8 yr" 1 may indicate an upper-limit to the sSFR of SMGs with a 
similar scaling to th e trends seen in lower activity galaxies at lower redshifts 
JPamen et alj2009h . 



burst phase. However, the dearth of SMGs at z < 1.5 with sSFR 
~ 10~ 8 -10~ 7 yr" 1 is not a selection effect and this upper enve- 
lope may be following the same trend in sSFR with r edshift seen 
in ga laxies with similar masses but lower SFRs (e.g iDamen et al.l 
120091) . 

Finally, we relate our new estimate of the redshift distribu- 
tion of SMGs to constraints on the e volution of their likely descen- 
dants: massive early- type galaxies JSwinbank etaUT 2006). As we 
have shown, the bulk of the SMG population with observed 870- 
/im fluxes above ~ 4mJy lie at redshifts of z ~ 1.5-3 with a 
median redshift of z ~ 2.5. We estimate that the volume density 
of SMGs at z = 2-3 above our flux is 1.2 x 10~ 5 Mpc" 3 , where 
we include both the identified and statistically identified samples in 
this estimate. Using a characteristic lifetime of the SMG phase of 
~ 100 Myr, we can correct this density for the burst duty cycle to 
derive a volume density for the remnants of 2 x 10" 4 Mpc " 3 . As 
we have shown, the estimated baryonic masses of these galaxies are 
~ 1.2 x 10 11 Mq combining ou r best estimate of the stellar mass 
with the typical gas masses from lGreve et al. I d2005h . If the burst of 
star formation we are seeing in the SMG phase is the last major star 
formation event in these galaxies then we expect their descendants 
to appear as passive, red galaxies at z ~ 1.5 (> 1 Gyr after z ~ 2). 

There have been various estimates of the volume density 
of massive, passive galaxies at z ~ 1-2 dMcCarthv et~al]|2004l ; 
iDaddi et alj2005tlTavlor et al.l2009al) . For galaxies with masses of 
> 10 11 Mq, the estimated space de nsities are 1-2 x 10" 4 Mpc -3 
(at z = 1.5-1.8;lTavlor et alj2009ah . 3 x 10" 4 Mpc" 3 (at < z >= 
1.7:lDaddi etdBpj) and 0.6 x 10" 4 Mpc ~ 3 (at < z >= 1.5; 
iMcCarthv et al.ll2004) . These estimates, with their various uncer- 
tainties, are comparable to the predicted volume density of massive, 
passive galaxies if these all undergo an SMG-phase at an earlier 
epoch. Hence, the starbursts in SMGs may be responsible for the 



formation of a large fraction of the passive, massive galaxies seen 
at 2 ~ 1.5. 

We can attempt a similar calculation comparing the SMG pop- 
ulation at 2 > 3 with the constraints on massive galaxies at z > 2. 
We estimate the volume density of z > 3 SMGs from i ]4.2l as 
2.8 x 10" 6 Mpc" 3 . This includes the 10 identified SMGs, 4 ± 2 
statistically identified SMGs and the remaining 21 ± 19 uniden- 
tified sources and assumes that they are contained within a red- 
shift range of z = 3-7. Using a characteristic lifetime of the 
SMG phase of ~ 100 Myr, we can correct this density for the 
burst duty cycle to derive a volume density for the remnants of 
3.8 x 10" 5 Mpc -3 . Again the estimated baryonic masses of these 
galaxies are ~ 1.2 x 10 11 Mq. Unfortunately observable limits on 
the volume density of passive galaxies are increasingly un certain 
at z > 2, but using the estimates from lCoppin et al] 120 091) of the 
volume density of massive galaxies of ~ 1-5 x 10" 5 Mpc" 3 , we 
again conclude that it is possible that the SMG population we have 
identified is also responsible for the formation of a significant of 
the most massive galaxies at z ~ 2.5. 

Thus we conclude that the presence of a sizable population 
of passive galaxies at high redshift may be intimately linked to the 
strong evolution in dust obscured starbursts in the distant Universe. 
Theoretical attempts to match the properties of high-redshift galax- 
ies therefore need to focus on these observable constraints as as- 
pects of the same problem dSwinbank et~aT]|2008f) . 



5 SUMMARY AND CONCLUSIONS 

We use deep multicolour imaging of the ECDFS in 17-bands to 
derive the photometric properties of the coun terparts of SMGs 
in the LESS LA BOCA survey of the ECDFS JWeifi et al.l [20091: 
Bigg s et al]|2010h . Our main results are as follows: 

(i) LESS radio, 24 and IRAC-identified SMGs have a me- 
dian redshift ofz = 2.2±0.1 and interquartile range of z = 1.8- 
2.7. Thus the peak activity in SMGs corresponds to the epoch of 
maximal quasar and star-formation activity in the Universe. The 
redshift distribution of LESS SMGs is con sistent with the spec- 
troscopic survey of [chapman et al. (2005), but higher than the 
photometric stu dies of the SHADES survey dOements et alj|2008l : 
iDve et al.l2 008). We find a higher-redshift tail to the distribution of 
LESS galaxies, with 10 (14%) identified SMGs at z > 3. Counter- 
parts identified through radio, 24 |j.m and IRAC emission have sta- 
tistically indistinguishable redshift distributions; similarly robust 
and tentative counterparts have comparable redshift distributions, 
albeit with some foreground contamination in the tentative sample. 
Previous studies provided tentative evidence that SMGs with the 
highe st submillimetre fluxes may be the highest redshift sources 
(e.g. Ilvison et alj|2002l ; iPope et al.ll2005l) . but with our extensive 
photometric redshifts we find no such correlation. 

(ii) A statistical study of the source population in the error cir- 
cles of the 55 SMGs that lack robust radio, 24 (xm and IRAC coun- 
terparts suggests that there is an excess of 26 ± 12 (2.2a) z > 1 
galaxies in these regions. Of these 14±8 (1.8<r) are atz — 2-3 and 
4 ± 2 (1.6cr) are at z — 4-5. This excess population corresponds to 
the counterparts or companions of the unidentified SMGs and our 
analysis then suggests that the redshift distribution of these uniden- 
tified SMGs peaks at z — 2.5 ± 0.3. This is similar to, but slightly 
higher than, that of the identified population, suggesting that many 
of the unidentified SMGs are at z ~ 2-3 and have radio or 24/^m 
fluxes just below our detection limits. 
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(iii) We estimate that there are 21 ± 19 (17 ± 15% of all the 
SMGs) LESS SMGs that are not robustly identified, and which are 
not accounted for in our statistical analysis of unidentified SMGs. 
These should be galaxies without any detectable mid-infrared emis- 
sion and as a result are likely to lie at z > 3. Including the identified 
SMGs at z > 3 we estimate that 28 ± 15% of all LESS SMGs lie 
at 2 > 3. 

(iv) We combine the redshift distribution of identified SMGs 
with that statistically determined for unidentified SMGs, including 
the SMGs which are not detected in our survey and are likely to lie 
at z > 3. We conclude that the likely median redshift of the entire 
population of 5870m'" ^ 4 mJy SMGs is z = 2.5 ± 0.6. 

(v) The separation of SMGs into those with redshifts above and 
below z = 1.4 broadl y agrees with the BzK colour-colour cri- 
teria jDaddi et al.ll20o4) . making the BzK colours a coarse but 
reliable redshift indicator for SMGs. Similarly the submillimetre- 
to-radi o flux ratios of LESS SMGs broadly agrees with the predic- 
tion of ICarilli & Y un (2000), although there is significant scatter, 
indicative of temperature variations between SMGs. This means 
that redshifts derived from submillimetre-to-radio flux ratios also 
exhibit significant scatter. Instead, we show that the IRAC colours 
of the ~ 90% of SMGs with z < 4 follow a trend with redshift of: 

2 = 2.1 + 1.91og 10 (S l 8/S , 3.6) and this may be a potentially useful 
redshift indicator with a \-a accuracy of a z = 0.4. 

(vi) The median rest-frame //-band absolute magnitude of the 
LESS SMGs is M H = -24.1 ±0.1 with an interquartile range 
of —24.7 to —23.6. Usi n g Mh and the average mass-to-light 
ratio from Hainline e t alj i 2010t) (converted to a Salpeter IMF) 
we calculate that the median stellar mass of the LESS SMGs is 
(9.2 ± 0.9) x 10 10 Mq, with an interquartile range of (4.7-14) x 
10 10 A/q. However, a x 2 analysis of the best-fit star-formation his- 
tories shows that, even with 17-band photometry spanning the ul- 
traviolet to mid-infrared, we cannot reliably distinguish different 
star-formation histories and ages for the SMGs. We estimate that 
this results in an additional factor of ~ 5 uncertainty in the mass- 
to-light ratios and hence the derived stellar masses. 

(vii) Using our photometric redshifts, submillimetre and radio 
fluxes we calculate that the median characteristic dust temperature 
of the SMGs is Td ~ 35.9 ± 1.4 K, with an interquartile range of 
28.5-43.3 K. The median far-infrared luminosity of the SMGs, de- 
rived from the radio luminosity, is Lfir = (8.2 ± 1.2) x 10 12 Lq 
and the interquartile range of Lfir = (3-13) x 1O 12 L . For a 
Salpeter IMF this corresponds to median SFR= 1100 MQyr -1 , 
with an interquartile range of 300-1900 MQyr -1 . We show that, 
for LESS SMGs, the apparent correlation between the far-infrared 
luminosity and Td is in part a selection effect. 

(viii) The far-infrared luminosity function of the LESS SMGs 
exhibits a strong redshift evolution, such that SMGs at z — 2- 

3 are more numerous and have higher luminosities than those at 
2 = 1-2. We find that the nor malisation of the lumin osity function 
is lower for LESS than for the lChapman et alJd2005h SMG sample, 
and by scaling the ECDFS submillimetre number counts we show 
that this is d ue to the underden sity of the ECDFS at submillimetre 
wavelengths JWeifl et all 20091) . 

(ix) The SFRD and fractional contribution to the global SFRD 
of the LESS SMGs with SsTO/am > 4mJy evolves with redshift 
and both peak at z ~ 2, where the LESS population contributes a 
total SFR density of 0.02 MQyr^Mpc -1 . If fainter submillime- 
tre sources have the same redshift distribution then SMGs with 
Ss70fim > 1 mJy produce ~ 50% of the SFRD of the Universe 
at 2 ~ 2. 

(x) The masses and the volume density of LESS SMGs at z = 



2-3 are comparable to those of massive, passive galaxies at z ~ 1- 
2, and similarly the volume density of z > 3 SMGs is comparable 
to the limits on the numbers of massive galaxies at z ~ 2-3. This 
suggests that a large fraction of the population of massive, passive 
galaxies at high-redshifts form most of their stars during an earlier 
SMG phase. 

This analysis demonstrates the strengths and weaknesses of 
photometric redshift analysis for SMGs in a field with excellent 
photometry. In the impending era of SCUBA-2 and Herschel 2> 
10° SMGs will be discovered; it will be impossible to obtain spec- 
troscopic redshifts for such large samples and hence the challenge 
is to obtain sufficient photometric coverage of these survey fields 
to allow a photometric analysis of the type described here. 
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Table 4. Observed photometry for robust counterparts to LESS SMGs. 3cr limiting magnitudes are presented where sources are covered by imaging but not detected; SMGs which are not covered by imaging in a 
given filter have no photometry listed in that filter. 
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LESS 2a 
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LESS 3 
LESS 6 
LESS 7 
LESS 9 
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LESS 10b 
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LESS 15 
LESS 16 
LESS 17 
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20.94 ± 
21.07 ± 

20.99 ± 
19.93 ± 
21.73 ± 

20.82 ± 

21.83 ± 
20.80 ± 
21.67 ± 
21.67 ± 

23.24 ± 
20.57 ± 
21.9^ 
21.2ffl% 
20.EftQt 
21.27^k 
20.6!J5t 
20.2»^t 
21.0©t 
21.8^t 
19.32^t 
20.73~jt 
20.8!>2t 
21.9$St 
21.81»}t 



0.09 
0.14 
0.09 
0.17 
0.05 
0.08 
0.09 
0.09 
0.08 
0.12 
0.08 
0.08 
0.05 
0.07 
0.06 
0.12 
0.09 
0.06 
0.08 
0.07 
0.17 
0.12 
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0.11 
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> 23.50 
21.33 ± 0.09 
19.67 + 0.04 
21.683E 0.10 
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APPENDIX A: DISCUSSION OF UNUSUAL SOURCES 

Most LESS SMGs have properties similar to previous SMG popu- 
lations. However, several robust counterparts have unusual proper- 
ties, X-ray emission, or 8 [im excesses indicative of the presence of 
an AGN. We discuss these galaxies on a case-by-case basis below. 
LESS 2a & 2b: LESS 2 has both a robust 24 ixm (LESS 2a) and a 
robust radio counterpart (LESS 2b), separated by 2.7", with z = 
1.80t°'i4 and 2 - 27 ±o.55 respectively. It is possible that the two 
counterparts are at the same redshift (z ~ 2), separated by ~ 20 
kpc and may in the process of merging. 

LESS 6: LESS 6 is ~ 1" from the robust 24ixm counterpart and 
~ 2.5" from the robust radio counterpart and has a photometric 
redshift of z = 0.40±g g§, therefore it is possible that the submil- 
limetre source is a background galaxy which is being gravitation- 
ally lensed by LESS 6. 

LESS 9: This SMG is X-ray luminous and has z phot = 4.63± ( j 1 ;^ 
making it the highest redshift X-ray source in our sample. Our best- 
fitting SED shows no 8 \im excess and thus no suggestion of an 



AGN from the near-infrared and optical photometry. 
LESS 10a: LESS 10a is one of two counterparts to LESS 10 iden- 
tified from extended or blended radio emission. It has some 8 \im 
excess over the best-fit SED, suggesting the presence of an AGN, 
but no detectable X-ray emission. 

LESS 11: Similarly to LESS 9 this galaxy shows no 8 ixm excess 
but is X-ray bright. 

LESS 19: Based upon the SED fitting LESS 19 has excess 5.8 
and 8 \im emission. However, the source is faint and the fit has 
X 2 = 10.0 so it is unclear whether the apparent excess is due to 
errors in the fitting or the presence of an AGN. 
LESS 20: This SMG is unusually radio bright with 5i.4ghz = 
4.25 mjy. There is no evidence of X-ray emission or an 8 ixm ex- 
cess, suggesting that LESS 20 is a radio-bright AGN. If this is the 
case then the AGN contribution to the radio flux means that Tb 
and Lfir are likely to be significantly over-estimated. Therefore, 
we exclude LESS 20 from our analyses of the SFRs and luminosity 
function of SMGs. 

LESS 40: LESS 40 has X-ray emission and an 8 \im excess and is 
highly likely to contain an AGN. 

LESS 50b: This galaxy is X-ray luminous but does not exhibit a 
compelling 8 ixm excess. 

LESS 57: LESS 57 has a strong 8 (im excess and is X-ray bright, 
compelling evidence for the presence of an AGN in this SMG. 
LESS 66: LESS 66 lies near the diffraction spike of a bright star, 
so some of the photometry may be unreliable. However, it is an 
optically bright point source with an 8 \un excess and X-ray emis- 
sion. There are broad emission lines in the spectra suggesting that 
LESS 66 is a submillimetre-bright quasar. 

LESS 67: This galaxy has coincident X-ray emission and may con- 
tain an AGN. 

LESS 74a & b: The two counterparts to LESS 74 have z = 
1.84±°;g and 1.71±g;?°, are separated by 2.7" (~ 20 kpc), and 
have some faint extended emission between them in the optical im- 
ages. Therefore, it is likely that LESS 74a and LESS 74b are under- 
going an interaction at z ~ 1.8 which triggered the submillimetre 
emission. 

LESS 75: Although LESS 75 is not X-ray detected there is strong 
8 and 5.8 ixm excess above the best-fit SED, indicative of a highly 
obscured AGN. 

LESS 84: LESS 84 is X-ray detected, and has a small 8 ixm excess 
suggesting it may contain an AGN. 

LESS 96: LESS 96 is similar to LESS 66 - it is X-ray luminous and 
has a strong 8 ixm excess above the best-fit SED, and it is an bright 
optical point source. We interpret this as evidence that LESS 96 is 
a submillimetre-bright quasar. 

LESS 106: This SMG is X-ray detected but no 8 ixm excess is ob- 
served. 

LESS 108: LESS 108 is identified as a bright local (z = 0.086) 
late-type spiral galaxy from its strong radio and 24 ixm emission. 
The similarity of the radio, mid- and near-infrared, and optical mor- 
phologies suggests that this is not a case of gravitational lensing. 
LESS 111: LESS 1 1 1 is X-ray detected and has a 8 \im excess, in- 
dicating it may contain an AGN. It lies < 3" from an extended 
foreground galaxy and is likely to be gravitationally lensed. 
LESS 114: This SMG is coincident with an X-ray source, but 
shows no evidence of an 8 [im excess above the best-fit SED. 
LESS 122: LESS 122 has excess flux at 8 ixm compared to the best- 
fit SED but is not X-ray detected. 
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APPENDIX B: SED FITS 

In Fig. lBll we show the measured photometry and best-fit SED for 
each SMG counterpart. The calculated photometric redshifts and 
errors are shown, and the probability distribution functions pre- 
sented for each galaxy. Most SMGs are well-fit by out template 
SEDs, although nine (12%) have excess 8 (im flux above the best- 
fit SED. 
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Figure Bl. Photometry and best-fit SEDs for robust SMG counterparts; points and error-bars show measured photometry and 
arrows represent 3cr detection limits. Error bars on the primary photometric redshift are 99% confidence limits. Inset panels show 
the minimum reduced x 2 at each redshift step with the photometric redshift primary and secondary solutions (where they exist) 
marked by solid and dashed lines respectively; the photometric redshift error is represented by the shaded region. 
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